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ABSTRACT 


In  this  thesis,  methods  for  forming  a  communications  link  between  Wireless 
Sensor  Networks  (WSNs)  by  enabling  each  WSN  to  act  as  a  smart  antenna  are  presented. 
Each  WSN  is  simulated  as  a  set  of  randomly  placed  sensor  nodes  within  a  planar  area. 
The  proposed  method  involves  a  searching  WSN,  a  receiving  WSN  and  a  link  budget  for 
establishing  the  link. 

The  searching  WSN  has  the  task  of  transmitting  a  search  beam  in  order  to  find 
adjacent  WSNs.  Like  a  lighthouse  this  is  done  in  a  rotating  beam  style  search  using  the 
sensor  nodes  as  an  aperiodic  array.  Results  show  that  for  a  random  array,  we  can  achieve 
a  specific  beamwidth  and  gain  as  a  function  of  the  number  of  elements  and  area.  We  also 
demonstrate  that  for  a  given  required  gain  level  we  can  spatially  thin  the  array  without 
significant  loss  of  gain  or  the  effects  of  grating  lobes. 

The  receiving  WSN  uses  a  spread  spectrum  based  space  division  multiple  access 
(SDMA)  receiver.  This  receiver  is  simulated  to  determine  the  direction  of  arrival  from  the 
searching  WSN  and  to  extract  the  location  information  from  the  searching  WSN’s  signal 
in  additive  white  Gaussian  noise.  From  the  DOA  and  the  location  information  within  the 
arriving  signal,  the  WSN  has  sufficient  knowledge  to  respond  to  the  query  of  the 
searching  WSN  and  form  the  communications  link. 
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EXECUTIVE  SUMMARY 


A  wireless  sensor  node  is  a  low  cost,  small,  battery  powered  electronic  device 
designed  to  monitor  or  measure  a  physical  phenomenon  of  the  environment  around  it, 
such  as  temperature,  pressure,  humidity,  soil  pH,  vibration,  motion,  light,  sound, 
radiation,  and  chemical  presence.  Wireless  sensors  are  generally  composed  of  a 
microprocessor,  the  appropriate  sensor(s),  a  transceiver  and  a  power  source.  Once  the 
sensor  makes  a  detection,  the  infonnation  collected  is  relayed  back  to  an  end  user  for 
further  processing  or  analysis.  A  collection  of  sensors  in  close  proximity  to  one  another 
may  form  a  Wireless  Sensor  Network  (WSN)  by  establishing  communications  with  one 
another  through  some  form  of  self-organization.  There  are  many  different  applications  for 
WSNs,  both  commercially  and  militarily.  In  this  work  we  assume  a  military  application 
where  the  sensor  nodes  are  deployed  in  hostile  territory,  and  therefore  delivery  of  the 
sensors  may  be  via  an  aerial  vehicle.  Delivery  in  this  manner  would  generally  lead  to  a 
random  spatial  distribution  of  the  sensors,  yet  the  sensor  nodes  are  close  enough  to  form  a 
local  wireless  sensor  network.  The  nodes  of  this  WSN  may  act  collectively  as  a  group 
for  such  actions  as  distributed  processing  or  forming  a  distributed  antenna  array. 

In  this  thesis,  we  assume  the  deployed  WSNs  must  transmit  their  information  to 
the  end  user  without  access  to  the  delivering  aerial  vehicle,  but  rather  terrestrially.  In  the 
solution  proposed  here,  where  the  WSNs  are  deployed  to  multiple  locations  nearby  to  one 
another,  they  can  fonn  a  terrestrial  Over  The  Horizon  (OTH)  network  by  forming  a 
multi-hop  communication  link  using  the  individual  WSNs.  In  order  to  form  this  link,  the 
WSNs  must  also  solve  the  additional  task  of  locating  one  another.  The  WSNs  act 
collectively  to  form  a  smart  antenna  array.  By  establishing  communication  links  between 
adjacent  WSNs,  an  extended  network  made  up  of  what  has  now  become  nodes  of  a  larger 
Over  The  Horizon  (OTH)  network  is  formed.  Once  the  OTH  network  is  connected,  the 
individual  WSNs  become  part  of  a  larger  network  that  can  relay  sensed  data  to  an  end 
user. 

The  goal  of  this  thesis  is  to  develop  and  validate  methods  to  establish  OTH 
communication  links  using  wireless  sensor  networks  by  enabling  each  to  act  as  a  smart 


antenna  array.  In  this  thesis,  we  approach  establishing  the  communication  link  between  a 
searching  WSN  and  a  receiving  WSN  as  a  four-step  process:  localization  of  the  sensor 
nodes  within  each  WSN  to  determine  its  geometry;  the  search  by  one  WSN  for  adjacent 
WSNs  by  acting  as  an  smart  antenna;  determination  of  the  direction  of  arrival  and  range 
of  the  searching  WSN  by  the  receiving  WSN;  and  then  forming  the  multi-hop  link  using 
WSNs  by  optimizing  the  gain  of  the  beam  forming  the  link. 

In  searching  for  the  adjacent  WSN,  each  of  the  sensors  act  in  concert  as  a  smart 
antenna  array.  To  describe  how  we  will  form  a  beam  using  a  group  of  randomly  placed 
sensor  nodes,  we  review  uniform  linear  and  planar  arrays  then  the  random  planar  array. 
Each  WSN  is  modeled  as  a  uniformly  random  distribution  of  sensors  nodes.  In  order  to 
link  the  WSNs,  we  will  wish  to  control  the  beamwidth,  gain  and  direction  of  the  array 
beam.  We  show  analytically  that  we  can  steer  the  beam  in  an  arbitrary  direction.  We  also 
demonstrated  control  over  the  beamwidth  and  gain  by  providing  MATLAB  results  and 
compare  the  results  to  analytical  solutions.  To  demonstrate  the  search  technique,  a  40 
element,  2.5k  x  2.5k  random  planar  array  scanning  the  horizon  is  presented. 

The  results  show  that  by  choosing  the  physical  size,  likely  as  sub-section  of  the 
WSN  we  can  choose  the  beamwidth  of  the  search  array.  Within  this  area,  we 
demonstrate  control  over  the  gain  by  choosing  the  number  of  elements  employed  in  the 
array.  Alternately,  for  a  given  required  gain  level  we  show  that  we  can  thin  the  array  (or 
reduce  density)  dramatically  without  significant  loss  of  gain  or  the  effects  of  grating 
lobes.  It  does  not  matter  which  particular  sensor  nodes  within  the  WSN  are  used  as  long 
as,  on  average,  they  are  randomly  distributed  within  the  physical  area  chosen.  This  is 
important  as  we  may  wish  to  distribute  the  burden  of  transmitting,  and  thus  consuming 
energy  among  the  nodes  within  the  WSN. 

From  the  perspective  of  the  receiving  WSN,  this  array  needs  to  know  the  direction 
of  arrival  (DO A)  and  range  to  the  searching  WSN  in  order  to  form  an  optimum  link.  A 
new  approach  to  DOA  determination  using  spread  spectrum  space  division  multiple 
access  (SDMA)  is  introduced,  demonstrated  for  a  uniform  linear  array  and  then  validated 
for  a  random  planar  array.  We  show  that  the  method  works  better  for  Walsh-Hadamard 


spreading  codes  than  for  pseudorandom  noise  codes  and  that  64  chips  is  sufficient  to 
accurately  establish  the  DOA  and  character  of  a  linear  antenna  array.  Using  a  random 
planar  array,  we  determined  the  DOA  of  an  incident  test  signal,  then  simulated  extracting 
a  message  data  stream  from  this  signal  in  a  noisy  channel  with  SNRs  of  20  db  and  -5  dB, 
accurately  extracting  the  message. 

A  link  budget  analysis  is  performed  for  a  simulation  involving  two  WSNs  using 
random  planar  arrays  on  a  flat  earth.  This  example  uses  commercially  available  sensor 
motes  at  two  transmit  power  levels.  Gain  calculations  are  made  for  each  power  level, 
solving  for  the  required  number  of  transmit  and  receive  elements  in  each  array.  In  this 
example  we  use  the  Crossbow  MICAz®  sensor  operating  at  2.4  GHz.  The  path  loss 
provided  the  greatest  amount  of  uncertainty  in  the  link  budget.  Of  the  models  examined, 
the  range  of  values  for  a  500  m  signal  at  2.4  GHz  varied  by  approximately  15  dB.  After 
determining  the  gain  of  the  random  planar  array,  we  found  that  for  a  500  meter  separation 
we  could  use  17  elements  on  the  receive  side  and  25  elements  on  the  transmit  side  when 
all  elements  were  operating  at  full  power.  When  all  were  operating  at  minimum  power, 
we  found  that  using  45  elements  to  receive  and  45  elements  to  transmit  provided 
sufficient  margin. 

In  this  thesis  we  have  developed  and  validated  methods  to  establish 
communication  links  between  wireless  sensor  networks  using  smart  antenna 
beamforming  techniques,  spread  spectrum  space  division  multiple  access  as  a  method  for 
a  receiving  WSN  to  determine  the  direction  of  arrival  and  message  sent  from  an  adjacent 
searching  WSN  and  the  use  of  link  analysis  to  optimize  the  number  of  sensors  used  to 
form  the  link. 
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I.  INTRODUCTION 


A  wireless  sensor  node  is  a  low  cost,  small,  battery  powered  electronic  device 
designed  to  monitor  or  measure  a  physical  phenomenon  of  the  environment  around  it, 
such  as  temperature,  pressure,  humidity,  soil  pH,  vibration,  motion,  light,  sound, 
radiation,  and  chemical  presence.  Wireless  sensors  are  generally  composed  of  a 
microprocessor,  the  appropriate  sensor(s),  a  transceiver  and  a  power  source.  Once  the 
sensor  makes  a  detection,  the  infonnation  collected  is  relayed  back  to  an  end  user  for 
further  processing  or  analysis.  Figure  la  shows  an  example  of  the  MICAz  sensor  mote  by 
Crossbow®  Technologies  [1]  and  in  Figure  lb  a  MICA2DOT  sensor  by  Crossbow  A 


(a)  MICAz  (b)  MICA2DOT  Sensor  Mote 

Figure  1.  Examples  of  Sensor  Motes  from  Crossbow  (From  [1]) 


A  collection  of  sensors  in  close  proximity  to  one  another  may  form  a  Wireless 
Sensor  Network  (WSN)  by  establishing  communications  with  one  another  through  some 
form  of  self  organization.  This  group  may  act  collectively  to  increase  sensing  or 
processing  power  by  performing  some  of  the  computations  onboard  the  sensor  node  itself 
[2]. 
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Many  different  types  of  applications  for  wireless  sensor  networks  exist,  in  both 
industry  and  the  military  world.  Because  they  are  wireless  and  self-organize,  wireless 
sensor  networks  have  the  potential  to  relay  information  without  a  large  infrastructure  cost 
or  physical  impact.  Commercial  applications  of  wireless  sensor  networks  include 
industrial  monitoring,  building  controls,  security,  traffic  management,  weather,  wildlife 
tracking,  and  agricultural  field  temperature-sensing  networks  [3],  Still  another  application 
consists  of  sensors  that  monitor  conditions  in  the  London  underground  tunnels  and  water 
systems  [4], 

Military  applications  may  include  remote  sensing  of  nuclear,  biological  or 
chemical  (NBC)  agents,  motion,  video,  etc.  It  is  assumed  that  the  sensor  nodes  are 
deployed  in  hostile  territory  and  therefore  delivery  of  the  sensors  may  be  via  an  aerial 
vehicle.  Delivery  in  this  manner  would  generally  lead  to  a  random  spatial  distribution  of 
the  sensors,  yet  the  sensor  nodes  are  close  enough  to  form  a  local  wireless  sensor 
network.  Once  established,  the  issue  becomes  how  to  get  the  information  from  the 
wireless  sensor  network  to  the  end  user.  The  nodes  of  this  network  may  act  collectively 
as  a  group  for  such  actions  as  distributed  processing  or  forming  a  distributed  antenna 
array.  Vincent  [5]  and  Chan  [6]  explored  the  use  of  beam  forming  to  communicate 
directly  with  an  overhead  UAV  as  a  method  to  relay  information  to  the  end  user. 

In  this  thesis,  we  assume  the  deployed  WSNs  must  communicate  their 
information  to  the  end  user  without  access  to  a  UAV.  In  the  solution  proposed  here, 
where  the  WSNs  are  deployed  to  multiple  locations  nearby  to  one  another,  they  form  a 
terrestrial  Over  The  Horizon  (OTH)  network  by  linking  the  individual  WSNs.  In  order  to 
form  this  link,  the  WSNs  must  also  solve  the  additional  task  of  locating  one  another. 
Each  WSN  acts  collectively  to  form  a  smart  antenna.  A  key  assumption  to  this  problem 
is  that  the  distance  between  WSNs  is  greater  than  the  individual  range  of  a  single  sensor 
node.  Thus  the  sensor  nodes  in  the  WSN  must  act  cooperatively  in  order  to  have 
sufficient  transmitting  gain  to  reach  the  next  WSN. 

This  problem  combines  and  draws  from  several  different  disciplines  in  order  to 
propose  a  solution.  Each  of  these  topics  will  be  discussed  further  in  this  thesis.  In  order 

to  form  the  initial  local  network,  the  concepts  of  localization  of  the  network  are  required. 
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Once  the  geometry  is  known,  the  sensor  nodes  of  the  WSN  act  collectively  to  form  a 
distributed  smart  antenna.  A  smart  antenna  generally  refers  to  an  array  antenna  that  can 
adapt  its  beam  pattern  using  a  computer  processor  [7].  Because  the  sensor  nodes  are 
scattered  in  a  random  manner,  the  concepts  of  random  processes  are  used  to  apply  them 
in  the  smart  antenna  composition.  The  concepts  of  spread  spectrum  and  signal 
processing  come  into  play  in  the  determination  of  the  direction  of  arrival  for  the  beams 
connecting  two  WSNs. 


A.  DESCRIPTION  OF  THE  PROBLEM 

We  begin  with  a  description  of  the  wireless  sensor  network,  depicted  in  Figure  2. 
Consider  a  set  of  wireless  sensor  nodes,  dispersed  within  an  arbitrary  boundary.  The 
sensor  nodes  form  an  ad  hoc  network  by  establishing  communications  with  one  another. 
The  boundary  is  defined  as  the  range  limit  of  the  communicating  sensor  nodes  within  the 
ad  hoc  network. 


Sensor 


Sensor  Node 
Connection 


WSN  Cluster 
Boundary 


Figure  2.  Wireless  Sensor  Network 


The  individual  sensor  nodes  may  be  randomly  positioned  within  a  boundary 
defining  the  WSN.  The  goal  of  this  thesis  is  to  explore  methods  for  forming 
communication  links  between  adjacent  WSNs.  By  establishing  communication  links 
between  adjacent  WSNs,  we  can  form  an  extended  network  made  up  of  what  has  now 
become  nodes  of  a  larger  Over  The  Horizon  (OTH)  network.  The  general  case  is  shown 
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in  Figure  3.  Once  the  OTH  network  is  connected,  the  individual  WSNs  become  part  of  a 
larger  network  that  can  relay  sensed  data  to  an  end  user. 


Sensed 

Phenomena 


WSN  D 
End  User 


WSN  A 


WSN  C 


WSN  B 


Figure  3.  Generalized  Over  The  Horizon  Network  of  Wireless  Sensor  Networks 


Next  we  consider  a  tactical  military  example  where  we  wish  to  link  WSNs  to 
form  an  OTH  network.  Consider  an  Unmanned  Aerial  Vehicle  (UAV)  or  other  aircraft 
that  disperses  sensor  networks  to  various  locations  in  hostile  territory,  as  envisioned  in 
Figure  4.  The  goal  is  to  connect  the  WSN  A  to  B,  B  to  C,  etc.  such  that  we  eventually 
have  a  link  from  A  to  E. 
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Figure  4.  Tactical  Application  of  Linked  Wireless  Sensor  Networks 

In  previous  work  related  to  this  scenario,  each  of  the  WSNs  communicated 
directly  to  the  UAV  when  it  was  available  overhead  [5],  In  this  work  we  seek  to 
communicate  not  via  a  UAV,  but  by  establishing  links  between  already  placed  WSNs. 

In  more  general  applications,  a  WSN  need  not  be  on  the  ground  but  can  be 
attached  to  a  building,  bridge,  structure  or  even  a  moving  vehicle.  The  goal  here  is  for 
the  WSNs  to  find  each  other  via  beam  forming  and  then  establish  a  communications  link 
at  the  highest  rate  possible. 


B.  THESIS  OBJECTIVE 

The  goal  of  this  thesis  is  to  develop  and  validate  methods  to  establish  OTH 
communication  links  using  wireless  sensor  networks  by  enabling  each  to  act  as  a  smart 
antenna  array. 
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To  begin  fonning  the  OTH  communication  link,  we  examine  methods  to  form  a 
search  beam  and  transmit  this  beam  around  the  horizon  in  order  to  illuminate  an  adjacent 
WSN.  We  will  assume  that  each  WSN  is  composed  of  a  uniformly  distributed  random 
array  of  sensor  nodes.  Using  theory  from  periodic  and  aperiodic  antenna  arrays,  we  will 
develop  methods  to  model  and  characterize  the  random  array  in  terms  of  gain  and  half 
power  beamwidth  as  a  function  of  the  number  of  elements  and  the  size  of  the  random 
array.  We  will  demonstrate  using  MATLAB  a  random  planar  array  forming  a  beam  and 
scanning  the  horizon. 

For  the  WSN  being  illuminated,  the  direction  of  arrival  (DO  A)  and  location  of  the 
searching  array  must  be  determined  in  order  to  form  an  optimum  link.  We  seek  to 
develop  and  simulate  methods  to  determine  the  DOA  and  location  information  using  a 
direct  sequence  spread  spectrum  based  space  division  multiple  access  technique  applied 
to  the  WSNs. 

In  order  for  the  WSN  to  act  as  a  smart  antenna,  the  local  geometry  of  each  WSN 
must  be  determined  in  order  to  define  the  sensor  node  positions  and  boundaries.  Because 
this  topic  has  been  extensively  reported  in  the  literature,  only  an  overview  of  applicable 
localization  techniques  is  presented  in  Appendix  A.  To  establish  an  optimum 
communications  link,  we  will  examine  a  pair  of  WSNs  using  a  link  budget  analysis.  We 
will  determine  the  gain  and  number  of  sensor  nodes  required  for  this  communication  link. 
Finally,  the  overall  OTH  network  is  realized  by  establishing  a  series  of  communication 
links  as  described  above.  This  area  is  present  conceptually  only,  as  the  focus  of  this  thesis 
is  to  investigate  the  establishment  of  the  link  between  two  adjacent  WSNs. 


C.  RELATED  WORK 

The  concept  of  beam  forming  using  an  array  of  randomly  distributed  sensor  nodes 
has  been  explored  previously  in  several  works.  Tong  [8]  studied  random  arrays  in  the 
context  of  randomly  distributed  micro-strip  patch  antennas  placed  on  the  sides  of  ships 
acting  as  a  digital  phased  array  radar  aperture.  Batson  [9]  studied  wireless  sensor 
networks  as  a  means  to  collect  signals  intelligence  by  building  adaptive  arrays.  Vincent 
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[5]  considered  a  similar  problem  to  the  one  we  explore  in  this  thesis,  namely  that  of  a 
WSN  deployed  in  a  remote  location  that  uses  adaptive  beamforming  to  communicate  to 
an  overhead  UAV.  Vincent  also  explored  energy  burden  sharing  across  the  sensor 
network,  a  concept  that  is  conceptually  useful  to  an  over  the  horizon  network.  Chan  [6] 
studied  adaptive  beamforming  for  direction  of  arrival  calculation  as  well  as  the  effects  of 
position  error  and  element  failure  on  the  beam  of  a  random  array  of  elements.  Godara 
[10]  presents  an  extensive  compilation  of  techniques  for  beamforming  and  direction  of 
arrival  estimation,  as  well  as  a  comparison  of  the  relative  performance  of  each.  While 
Hong,  et  al  [1 1]  use  cooperative  methods  to  establish  a  relay  network  of  sensors,  they  use 
spatial  diversity  techniques  rather  than  beamforming. 

This  thesis  builds  upon  the  concept  of  the  random  planar  array  of  sensor  nodes, 
but  here  the  objective  is  to  build  an  over  the  horizon  network  of  WSNs.  We  also  propose 
new  ways  for  searching,  tracking  and  establishing  a  link  between  adjacent  WSNs.  The 
communication  links  between  adjacent  WSNs  are  established  by  using  beam  searching 
and  direction  of  arrival  methods  based  upon  spread  spectrum  techniques. 


D.  THESIS  ORGANIZATION 

This  thesis  is  organized  as  follows.  Chapter  II  provides  a  link  budget  analysis  to 
determine  the  gain  requirements  for  the  transmitter  and  receiver  given  data  rates,  range 
and  power  considerations.  Chapter  III  describes  the  use  of  beam-forming  as  a  method  to 
search  for  adjacent  WSNs.  Simulation  results  for  random  planar  arrays  and  search 
patterns  are  included.  Chapter  IV  describes  the  link  establishment  process  and  the 
determination  of  the  direction  of  arrival  (DOA)  using  spread  spectrum  based  space 
division  multiple  access  (SDMA)  methods.  Simulation  results  of  a  random  planar  array 
detecting  the  arrival  of  a  signal  of  interest  as  well  as  extracting  a  sample  data  stream  are 
presented.  Chapter  V  presents  the  conclusions  of  this  work  and  topics  for  further 
research.  Finally,  Appendix  A  summarizes  the  sensor  node  localization  methods  and 
Appendix  B  includes  MATLAB  code  used  in  this  thesis. 
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II.  COMMUNICATIONS  LINK:  LINK  BUDGET 
CONSIDERATION 


The  over  the  horizon  (OTH)  communications  system  is  essentially  a  network  of 
networks  linking  the  individual  WSNs.  Several  issues  must  be  resolved  in  order  to  link 
these  WSNs  of  an  OTH  network.  This  thesis  examines  the  link  between  just  two  of  the 
WSNs,  but  for  completeness,  we  discuss  the  bigger  issues  of  the  OTH  network  as  a 
whole.  Many  of  these  issues  are  not  too  different  from  the  issues  discussed  in  Appendix 
A  on  localization.  In  this  chapter,  we  begin  with  a  conceptual  look  at  the  OTH  network, 
then  consider  just  two  WSNs  of  the  overall  chain  of  WSNs  that  eventually  form  the  OTH 
network.  Within  this  view,  we  consider  them  as  two  points  of  a  communications  link  and 
apply  a  link  budget  analysis. 


A.  BUILDING  THE  OVER  THE  HORIZON  NETWORK 

Each  WSN  of  the  OTH  network  must  find  the  adjacent  WSN  to  connect  to.  Once 
the  initial  link  is  established,  say  from  A  to  B  as  in  Figure  3,  WSN  B  must  now  continue 
the  link  and  connect  to  a  third  WSN,  for  example  WSN  C  of  Figure  3.  While 
establishing  the  second  link  the  WSN  node  must  still  maintain  contact  with  the  first 
WSN.  A  particular  WSN  may  subdivide  as  a  method  to  establish  contact  with  more  than 
one  other  WSN.  The  problem  of  subdividing  an  individual  WSN  may  be  addressed  in 
several  ways.  Perhaps  the  WSN  could  physically  partition  the  sensor  nodes  to  address 
and  connect  to  other  nodes,  e.g.  the  Eastern  most  nodes  form  an  array  to  communicate  to 
another  WSN  to  the  East,  and  the  Western  most  sensor  nodes  form  an  array  to 
communicate  to  the  West,  as  depicted  in  Figure  5.  In  Figure  5,  we  show  WSN  B 
subdivided  laterally,  such  that  the  red  nodes  form  an  array  and  communicate  with  WSN 
A  and  the  green  nodes  form  an  array  and  communicate  with  WSN  C. 
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Figure  5.  WSN  B  Sensor  Nodes  Partitioned  to  Simultaneously  Communicate 

With  WSN  A  and  WSN  B 


Alternately,  the  WSN  could  employ  some  subset  of  the  sensors  using  a 
multiplexing  technique  based  on  time,  code  or  space.  Additionally,  the  frequency  band 
could  be  split  into  different  frequency  channels  for  communicating  to  individual 
surrounding  WSNs.  This  thesis  does  not  address  the  details  of  subdividing  the  WSNs,  but 
rather  approaches  this  part  of  the  overall  process  conceptually. 

In  the  following  chapters,  we  will  model  the  transmitting  and  receiving  WSNs  as 
smart  antennas,  taking  advantage  of  the  element  localization  information,  the  random 
nature  of  the  element  spatial  distribution  to  form  a  transmit  beam  and  a  spread  spectrum 
SDMA  technique  to  determine  the  relative  geometries  between  WSNs. 

Essentially,  the  WSN  will  now  be  functioning  as  individual  antennas,  and  thus  the 
communication  link  can  broadly  be  described  by  Figure  6.  This  figure  depicts  a 
terrestrial  communications  link  formed  over  some  distance  R  between  antennas  A  and  B. 
Because  of  the  nature  of  the  WSNs,  we  can  change  some  of  the  characteristics  of  the 
communication  link,  such  as  its  effective  aperture  and  the  number  of  elements  used  to 
form  the  antenna.  This  chapter  seeks  to  examine  those  parameters  given  the  known 
parameters  such  as  range,  required  data  rate,  path  loss  and  available  transmitter  power. 
We  begin  with  a  discussion  of  the  link  budget  equation  and  then  perform  an  analysis 
using  a  commercial  sensor  mote  and  an  example  random  planar  array. 
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Figure  6.  Simplified  Terrestrial  Radio  Communications  Link 


B.  LINK  BUDGET  EQUATION 

The  link  budget  allows  us  to  perform  a  tradeoff  analysis  when  we  are  given 
control  over  components  of  the  link.  In  this  section,  we  examine  the  individual  terms  of 
the  link  budget  equation  [12]: 


(1) 


V  o  J  Re  qd 


where  Mcib  is  the  link  margin  (dB),  P,  is  the  transmitter  power  (dBW),  G,  is  the  gain  of  the 
transmitting  array  (dBi),  G,  is  the  receiving  array  gain  (dBi),  Eh  /  N0  is  the  energy  per  bit 

to  noise  power  spectral  density  ratio  (dBW/Hz),  R/>  is  the  bit  rate  of  the  communications 
between  WSNs  (dB-bit/sec),  k  is  Boltzman’s  constant  (1.38  x  10  23  J/K),  Ts  is  the 
effective  system  temperature  (K°)  generally  given  as  Ts  =  Ta  +  Te , where  Ta  is  the  antenna 

noise  temperature  and  Te  is  the  receiver  noise  temperature.  The  path  loss,  Ls  (dB)  is  the 
loss  of  a  signal  over  the  distance  from  point  A  to  point  B,  including  multipath,  fading, 
two  ray  bounce,  weather,  etc.  and  L0  includes  “other  losses”  such  as  transmitter 
inefficiencies,  line  loss,  polarization  mismatch,  etc.  [12].  The  link  margin  is  the 
difference  between  what  we  need  to  establish  the  link  and  what  we  have  available  from 
the  combined  components  of  the  communications  link.  We  assume  that  each  transmitter 
can  vary  its  power  output  in  discrete  steps,  but  that  all  elements  transmit  at  equal  levels. 
The  gain  of  the  transmitting  and  receiving  arrays  refer  to  the  gain  of  the  random  planar 
array  of  judiciously  selected  elements.  The  gain  of  a  random  array  of  transmitters  will  be 
described  in  Chapter  III  on  random  arrays  and  will  be  one  of  the  key  parameters  we  can 
vary.  Here,  we  will  assume  the  “other  losses”  are  negligible. 
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1.  Link  Analysis 

Here  we  consider  two  WSNs  deployed  into  separate  but  nearby  regions  on  a  flat 
environment,  not  unlike  the  scenario  depicted  in  Figure  4.  For  example,  let  us  consider 
as  sensors  the  Crossbow'  MICAz,  a  commercially  available  sensor  mote,  to  apply  to  this 
scenario.  MICAz  operates  at  2.4  GHz  and  generates  RF  transmitter  power  from  40 uW  to 
hnW.  The  nominal  outdoor  range  for  a  single  receiver/transmitter  pair  is  given  as  100m 
at  the  receive  sensitivity  of -120  dBW  [13].  The  maximum  data  rate  for  this  mote  is  250 
kbps.  Using  devices  like  these  one  may  create  a  WSN  containing  50  sensors  arrayed  over 
a  5  m  x  5  m  (40/1  x  40/1)  region.  For  simplicity,  we  will  assume  that  both  WSNs  are 
the  same  size  and  contain  the  same  number  of  elements  and  that  the  separation  distance 
between  the  two  WSNs  is  500  m  (4000  A). 

We  will  assume  that  all  50  of  the  receive  sensor  nodes  are  awake  and  in  “receive” 
mode,  although  generally  this  need  not  be  the  case.  We  now  evaluate  the  terms  in  the 
link  budget  equation  for  this  example.  The  digital  modulation  scheme  for  the  MICAz  is 
BPSK  or  OQPSK  per  the  IEEE  803.15.4  standard.  As  a  first  example,  we  assume  that 
each  transmitter  is  operating  at  the  maximum  power  level  of  1  rnW  which  equates  to  -30 
dBW.  At  minimum  power,  the  transmitters  produce  40  juW  or  -54  dBW  which  we 
consider  second.  We  assume  that  all  the  transmitters  operate  at  the  same  power  level.  For 
the  path  loss  term  of  the  link  budget,  the  contribution  is  not  as  clear.  Before  continuing 
with  the  link  equation,  we  briefly  examine  the  path  loss  for  two  sensor  nodes  positioned 
on  the  ground. 


a.  Path  Loss  Models 

Many  models  exist  for  land  mobile  radios  and  wireless  LAN  (WLAN) 
systems  and  produce  a  wide  range  of  values  when  applied  to  sensor  nodes.  A  model  to 
estimate  path  loss  for  sensor  nodes  at  ground  level  is  not  readily  available  as  most 
address  antenna  heights  of  approximately  1  m  or  above.  Models  that  do  include  antenna 
height  are  relatively  sensitive  to  the  height  of  the  transmit  and  receive  antennas.  For  the 
MICAz  sensor  mote,  which  uses  a  one-half  wavelength  dipole  (~6.25  cm)  and  sits  on  a 
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battery  pack  and  circuit  board  (~3  cm),  this  gives  us  an  approximate  antenna  height  of  10 
cm.  We  now  examine  the  most  common  models  for  path  loss. 

The  free  space  path  loss  at  a  distance  d  and  wavelength  X  is  given  by  [12]: 

t  im  r 4tt d\2 
4  =  101Og,o  -f” 

V  A 

The  free  space  path  loss  is  not  applicable  to  the  scenario  under 
consideration,  but  it  serves  as  a  baseline  for  comparison  with  the  other  models  that  will 
be  presented  below.  In  this  case,  the  antennas  are  very  close  to  the  ground  and  therefore 
the  path  loss  is  expected  to  be  greater.  Green  and  Obaidat  developed  a  path  loss  model 
for  wireless  LANs  operating  at  2.4  GHz  that  takes  antenna  height  into  account  and  was 
experimentally  studied  under  conditions  similar  to  ours.  The  Green  and  Obaidat  path  loss 
equation  is  given  by  [14]: 

lwlan  =  7.6  +  40  log10  d  -  20  log10  h,hr  ( dB ) 

where  d  is  the  separation  distance  and  h,  and  hr  are  the  heights  of  the  transmit  and  receive 
antennas,  respectively. 

Okumura  [15]  and  Hata  [16]  produced  landmark  papers  on  propagation 
loss  in  land-mobile  radios,  but  these  models  are  for  signals  below  1500  MHz  and 
generally  for  antenna  heights  of  1  m  or  more.  The  Green  and  Obaidat  models 
demonstrated  very  good  agreement  with  their  experimental  results,  for  antenna  heights  of 
1-2  m,  but  produce  higher  losses  when  used  for  antenna  heights  of  5-10  cm. 

A  propagation  model  based  on  a  flat  Earth  assumption,  where  the  Earth  is 
represented  as  a  dielectric  is  another  approach.  Linnartz  [17]  provides  the  path  loss  at  a 
distance  d  for  the  flat  Earth  as: 

r  A2  r  .  (IK  h  VI 

LnF  =  - - -y  2  sin  ~—hth 

pE  {4Kd)2[  \Xd  r)\ 

where  X  is  the  wavelength,  and  ht  and  hr  are  the  heights  of  the  transmit  and  receive 
antennas,  respectively.  Because  the  antennas  are  so  close  to  the  ground  and  we  are  using 
a  2.4  GHz  example  this  model  would  only  apply  in  a  very  limited  number  of  cases  where 
the  environment  was  free  of  scattering  objects. 
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The  calculated  path  losses  for  the  2.4  GHz  signal  at  a  separation  of  500  m 
with  both  antennas  at  10  cm  are  listed  as  follows: 


Free  Space 

94  dB 

Hata  Model 

121  dB 

Green/Obaidat  Model 

134  dB 

Based  upon  the  operating  environment,  the  path  loss  for  this  system  ranges  from  94  to 
134  dB.  Each  of  the  models,  except  Free  Space,  is  sensitive  to  antenna  height.  The  Green 
and  Obaidat  model  seems  reasonable  for  use  in  this  work  given  that  their  model  was 
developed  specifically  for  WLAN  applications  at  2.4  GHz  and  had  good  agreement  with 
experimental  results.  However,  this  model  used  much  greater  antenna  heights  suggesting 
that  using  this  model  for  sensor  nodes  on  the  ground  may  not  produce  accurate  results. 

b.  Link  Budget  Gain  Calculations 

We  wish  to  know  what  level  of  gain  is  required  for  the  transmit  and 
receive  antenna  arrays  forming  the  communications  link  shown  in  Figure  6.  We  can  solve 
the  link  equation  to  determine,  with  margin  M,  what  gain  the  two  WSN  arrays  must 
produce.  Since  we  are  given  the  receive  power  sensitivity  Pr,  of  the  MICAz,  assumed  to 
incorporate  the  data  rate  and  required  Eb/Na  ,  Equation  (1)  simplifies  to: 

MdB  =  Pt+Gt  +  Gr—Pr  —  Ls  (2) 

where  Mjb  represents  the  difference  between  the  actual  received  power  and  the  minimum 
useable  power,  given  by  the  receive  sensitivity.  Also  note  that  the  total  power 
transmitted  by  the  array  of  Nt  transmitting  elements  is  N,Pt . 

Using  Equation  (2)  and  combining  the  tenns  from  the  link  analysis  and  the 
Green/ Obaidat  path  loss  estimate,  we  have  for  the  1  mW  (-30dBW)  transmitter: 

M  =  (Gt  +  Gr)  +  Nt  +Pt-Pr-Ls 
M  =  (Gt  +  Gr)  +  N,  -  30  -  (-120) -  134) 

(G,  +  Gr)  =  M-N,  +  44 

Thus  the  combined  gains  of  the  transmitting  array  and  receiving  arrays 
must  meet  a  margin  balanced  by  the  total  transmitted  power  less  a  factor  for  path  loss. 


14 


Excluding  the  free  space  model,  the  range  of  path  loss  values  is  on  the  order  of  15  dB, 
suggesting  that  we  should  choose  this  as  our  margin. 

When  the  transmitters  are  operating  at  minimum  power  (-54  dBW),  from 
Equation  (2)  the  required  gain  is: 

(G,+Gr)  =  M-N,+  68 

In  this  case,  it  is  clear  that  we  must  now  either  decrease  our  acceptable  margin  or  increase 
the  number  of  transmitters  or  both  in  order  to  form  the  communications  link.  Because 
transmission  consumes  much  more  energy  than  reception,  it  is  envisioned  that  we  will 
use  many  more  receive  nodes  than  transmit  nodes  to  meet  the  gain  requirement.  For 
example  we  may  use  15  transmit  and  30  receive  nodes.  In  Chapter  III,  we  will  be  able  to 
determine  these  numbers. 

Clearly  the  path  loss  dominates  the  overall  margin  requirement.  One 
method  to  reduce  the  path  loss  is  to  elevate  the  sensor  nodes.  This  increases  the  antenna 
height  and  rapidly  reduces  the  path  loss,  thus  leading  to  greater  range.  Typically, 
doubling  the  antenna  height  increases  the  received  power  by  a  factor  of  4  [17].  For 
example,  the  sensor  network  could  be  placed  on  the  side  of  a  building  or  rooftop, 
elevated  on  top  of  a  vehicle,  attached  to  lamp  posts,  etc. 

In  summary,  we  applied  the  link  budget  equation  to  a  link  formed  between 
two  WSNs  and  determined  the  required  combined  gain  from  the  receive  and  transmit 
antennas.  In  the  next  chapter  we  examine  how  to  calculate  this  required  gain  as  well  as 
some  other  parameters  of  the  antenna  we  will  form  from  the  WSN. 
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III.  THE  SEARCH  FOR  OTHER  WIRELESS  SENSOR 

NETWORKS 


In  this  chapter,  we  approach  the  communications  link  from  the  viewpoint  of  a 
searching  WSN  in  the  eventual  OTH  network.  For  example,  this  could  be  WSN  A  of 
Figure  3,  which  has  sensed  a  phenomenon  and  seeks  to  pass  this  information  along  to  the 
end  user.  This  WSN  has  the  task  of  transmitting  a  search  beam  in  order  to  find  adjacent 
WSNs.  Like  the  rotating  beam  of  a  lighthouse,  the  WSN  will  steer  its  beam  in  a  rotating 
manner  around  the  horizon  by  using  the  local  network  of  sensor  nodes  as  an  aperiodic 
antenna  array. 

In  order  to  address  the  method,  we  first  discuss  the  concepts  of  periodic  arrays. 
Then  we  extend  these  concepts  to  the  random  array,  a  type  of  aperiodic  array.  Next  we 
explore  the  characteristics  of  the  two  dimensional  random  array,  using  MATLAB  to 
calculate  gain,  beamwidth  and  sidelobe  levels.  Finally,  we  discuss  the  issues  of  the 
search  between  WSNs,  using  random  planar  arrays. 


A.  BEAM  FORMING  WITH  PLANAR  ARRAYS  OF  RANDOMLY  PLACED 

SENSOR  NODES 

To  describe  how  we  will  form  a  beam  using  a  group  of  randomly  placed  sensor 
nodes,  we  begin  by  reviewing  the  theory  of  the  uniform  linear  array  (ULA)  and  the 
uniform  planar  array.  From  this,  we  make  the  extension  to  planar  arrays  containing 
randomly  placed  elements.  In  each  case,  we  assume  that  the  sensor  antenna  elements  are 
half  wave  dipoles  oriented  vertically  to  the  plane. 

1.  Geometric  Arrays 

The  far  field  pattern  for  an  array  of  geometrically  arranged  antennas  is  calculated 
using  pattern  multiplication;  that  is,  the  far  field  radiation  pattern  F(0)  can  be  expressed 
as  the  product  of  the  element  pattern  E0(6)  and  the  array  factor  AF{0)  [18]: 

F(8)  =  Ea(6)xAF(e) 
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For  a  z  -directed,  half  wavelength  dipole,  the  normalized  element  factor  is  [18]: 


cos 


Km- 


K  ' 
cos  0 

2 


sin  6 


(3) 


and  is  assumed  to  be  the  element  in  all  the  gain  pattern  calculations  that  follow.  The 
array  factor  characterizes  the  radiation  pattern  based  upon  the  geometry  of  the  antenna 
elements.  We  now  discuss  the  one-dimensional  uniform  linear  array. 


a.  Uniform  Linear  Array  (ULA) 

Figure  7  shows  the  geometry  of  a  uniform  linear  array  and  its  coordinate 
system,  and  observation  direction  R(0).  The  angle  0  is  measured  from  the  z  axis. 


Z 


Figure  7.  Uniform  Linear  Array  Geometry  and  Coordinate  System 

For  an  equally  spaced  set  of  N  isotropic  radiators  with  wave  number 
k-2nlX,  spaced  by  a  distance  d  along  the  x  axis,  the  array  factor  may  be  expressed  as: 
[18] 

N-l 

AF(0)  =  ^Anejnkdsiae  (4) 

n= 0 
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where  An  is  the  amplitude  of  the  drive  of  the  /2th  current  element,  and  is  the  phase 

contribution  of  nth  element.  Note  that  the  array  factor  is  rotationally  symmetric  around 
the  x-axis.  Figure  8  illustrates  how  the  observation  angle  6  and  the  spacing  d  affect  the 
phase  term  of  the  array  factor.  For  example,  when  6  is  zero,  all  the  elements  have  the 
same  phase  because  sin($)  is  zero;  when  sin(6>)  =  1,  the  phase  contribution  from  each 
element  is  equal  to  kd.  From  this,  we  see  that  the  spacing,  observation  angle  and  wave 
number  determine  the  phase  of  the  far  field  array  factor. 


Element 

Figure  8.  Phase  Difference  is  determined  by  0  and  spacing  d 

We  can  simplify  Equation  (4)  by  letting 

i//  =  kd  sin  0 

and  letting  all  the  current  amplitudes  be  equal,  i.e., 

4  =  4  =  -  =  4-i 

Thus,  Equation  (4)  becomes: 

AF(y,)  =  a/X>  (5) 

n= 0 

For  N  elements,  and  the  beam  pointing  in  the  0=0  direction,  the  array  factor  reduces  to 

AF{¥)  =  A0N 

Thus,  the  maximum  array  factor  is  proportional  to  N. 
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It  can  be  shown  [18]  that  Equation  (5)  can  be  expressed  as  the  more 
common  expression  for  the  normalized  array  factor  / ,  or  pattern  factor  of  a  uniformly 


excited,  equally  spaced  linear  array,  centered  at  the  origin: 

sin(Ay/2) 

v  '  Nsm(yf/2) 


(6) 


A  close  examination  of  Equation  (6)  reveals  some  important  characteristics  of  array 
antennas  [18]: 


•  As  the  length  of  the  array  ( Nkd)  increases,  the  width  of  the  main  lobe 
decreases,  i.e.,  the  half  power  beamwidth  (HPBW)  becomes  smaller. 


•  As  N  increases,  the  side  lobe  levels  decrease.  The  Side  Lobe  Level  (SLL) 
is  expressed  as  the  ratio  of  the  maximum  value  of  the  largest  side  lobe  to 
the  maximum  value  of  the  main  lobe.  This  is  usually  described  in  the 
literature  as  “dB  below  the  main  lobe.” 


A  characteristic  of  periodic  arrays  is  revealed  when  the  spacing  between 
elements  increases  beyond  a  certain  level,  and  energy  from  the  array  is  directed  into 
unintended  directions.  These  unintended  beams  from  the  array  essentially  reduce  the  gain 
of  the  array  along  the  desired  direction  by  putting  energy  into  unwanted  directions. 
These  unintended  beams  are  called  grating  lobes  and  are  the  subject  of  the  next  section. 


b.  Grating  Lobes 

Grating  lobes  occur  when  element  spacing  is  greater  than  one-half 
wavelength  [18].  These  lobes  are  often  of  equal  intensity  to  the  main  lobe(s)  of  the  array 
and  are  due  to  the  periodicity  of  the  element  spacing.  Grating  lobes  do  not  occur  for 
arrays  with  spacing  less  than  one-half  wavelength;  thus,  nearly  all  periodic  arrays 
conform  to  this  spacing  constraint  [18].  In  order  to  demonstrate  the  occurrence  of  grating 
lobes,  Figure  9  shows  the  pattern  factor  for  an  eight  element  ULA,  with  a  spacing  of  0.5, 
1.0  and  1.25  wavelengths  between  elements.  Only  the  first  quadrant  is  shown,  as  the 
pattern  is  symmetric  about  0°.  The  blue  line  shows  the  pattern  factor  for  an  array  with 
one-half  wavelength  spacing  and  no  grating  lobes.  The  red  and  green  lines  show  the 


20 


pattern  factor  for  arrays  with  a  spacing  of  1.0  and  1.25  wavelengths  and  the  resultant 
grating  lobes.  Notice  that  for  the  one-half  wavelength  spaced  array,  shown  in  blue,  most 
of  the  energy  of  the  array  is  concentrated  in  the  main  beam,  while  for  the  larger  spacing 
the  energy  is  dispersed  into  the  grating  lobes. 


Figure  9.  Pattern  Factor  for  an  8  Element  ULA  Showing  Grating  Lobes  as  a  Result  of 

Spacing  Greater  Than  One-Half  Wavelength 

Thus  the  designer  is  presented  with  a  dilemma  -  in  order  to  have  a  small 
beamwidth,  the  ULA  must  be  physically  large,  but  increasing  the  length  will  necessarily 
result  in  either  a  higher  complexity  due  to  the  greater  number  of  elements  or,  if  spread 
out  too  much  to  keep  the  number  of  elements  down,  grating  lobes.  It  will  be  shown  later 
that  since  random  arrays  do  not  have  a  periodicity,  grating  lobes  generally  do  not  occur. 

We  now  extend  the  concepts  of  the  uniform  linear  array  to  form  a  two- 
dimensional  periodic  array. 

c.  Two-Dimensional  Periodic  Arrays 

Figure  10  shows  the  geometry  for  a  two-dimensional  periodic  array. 
Elements  are  geometrically  arranged  in  the  x-y  plane,  and  the  far  field  observation 
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direction  is  shown  by  /?(0,0).  For  this  array,  the  principle  of  pattern  multiplication  is 
also  applied.  Here  we  can  write  the  far  field  pattern  F(Q,(p)  =  Fc(0,0)  x  Fv(0,0) ,  where 
Fx(6,(j))  is  the  pattern  of  a  linear  array  along  the  x  axis,  and  FJO.cf))  is  the  pattern  of  a 
linear  array  along  they  axis  [19]. 


Z 


Figure  10.  Uniform  Two  Dimensional  Array  Geometry  and  Coordinate  System 


The  array  factor  for  a  two  dimensional  periodic  array  with  M  elements  of 
spacing  dx  in  the  x  direction  and  N  elements  of  spacing  dv  in  the  y  direction  is  given  by 
[20]: 


AF(0,0)  =  ££a„ 


gjk(m—\)(dx  sinflcos tp+fix  )^Mn-l)(dy  sin0sin0+j3v  ) 


m= 1 n= 1 


(7) 


Combining  the  exponential  terms  yields: 


M  N 


AF(e,^  =  JJJJA,me 


jk[(m-l)(dx  sin0cos0+/3x )+(«— l)(dy  smdsm(j)+py)] 


(8) 


where  Amn  is  the  current  distribution  at  the  (m,n)th  element.  The  array  can  be  steered  to  an 
arbitrary  direction  (0o,(j)o)  via  setting  [20]: 
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f5x  =  -kd  sin 0O  cos  0O  and  /3,  =  -kd  sin0o  sin 0o  (9) 

When  combined  with  a  specific  radiation  element,  the  total  radiation 
pattern  is  found  by  multiplying  the  element  pattern  Fo(0,0)  with  the  two  dimensional 
array  factor  AF(0,0)  [18]: 

F(0,0)  =  Fo(0,0)xAF(0,0)  (10) 

The  element  pattern  is  the  radiation  pattern  of  a  single  element  in  0  and  (f >.  In  our  case,  the 
radiation  pattern  for  a  half  wavelength  dipole  given  is  rotationally  symmetric  in  the  angle 
0  ;  therefore,  Equation  (3)  can  be  used  in  Equation  (10)  as  E0  (0,0) . 

It  can  be  seen  that  with  half-wavelength  spacing,  even  a  modest  7x7 
wavelength  array  would  contain  nearly  50  elements.  It  is  possible  to  thin  out  these  arrays 
yet  still  have  control  of  the  beamwidth,  gain  and  grating  lobes.  One  can  thin  the  arrays 
by  placing  the  elements  randomly  within  a  given  area  at  a  much  lower  density  than  that 
of  a  periodic  array.  The  next  section  develops  these  concepts. 

2.  Random  Arrays 

A  random  array  is  an  array  antenna  whose  elements  are  no  longer  defined  by  a 
fixed  geometric  spacing  but  rather  the  element  locations  are  now  random  variables. 
Figure  1 1  depicts  the  periodic  linear  array,  the  random  linear  array  and  the  random  planar 
array.  A  one-dimensional  periodic  linear  array  has  a  fixed  inter-element  spacing,  while 
the  element  spacing  in  a  random  linear  array  is  a  random  variable.  The  elements  are 
randomly  placed  along  both  the  x  and  y-axes  in  a  random  planar  array,  i.e.,  the  inter¬ 
element  spacing  is  random  along  both  axes. 
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Periodic  Linear  Array 


Random  Linear  Array 
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Random  Planar  Array 


Figure  1 1 .  Illustration  of  Periodic  Linear  Array,  Random  Linear  Array  and  Random 

Planar  Array  Showing  Example  Element  Distribution 

The  array  factor  for  a  random  linear  array  is  similar  to  the  array  factor  for  a 
periodic  array,  repeated  here  for  convenience: 

N- 1 

AF0(Q)  =  £  Aneikndsme 

n= 0 

except  that  now  the  element  spacing  term  in  the  exponential  is  a  random  variable.  For 
the  purposes  of  performance  comparison,  the  periodic  array  is  termed  the  design  array 
and  its  array  factor  is  the  design  array  factor  AF0(6 ) .  The  phase  contribution  from  each 
element  in  a  random  linear  array  is  no  longer  a  deterministic  term,  i.e.,  ndx ,  but  rather 
from  the  actual  x  coordinate  of  the  nth  element,  xn,  and  thus  the  phase  is  also  a  random 
variable.  The  array  factor  for  linear  random  array  becomes: 

jV— 1 

AF(Q)  =  ^  AnejkXn  sine  (11) 

n= 0 

Because  the  element  location  is  now  a  random  variable  the  array  factor  is  also  a 
random  variable.  It  is  customary  to  represent  the  array  factor  of  a  random  array  in  terms 
of  the  ensemble  average,  referred  to  as  the  average  array  factor,  given  by: 

AF(d)  =  E[AF(6)] 
where  £[•]  is  the  expectation  operator. 
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As  two  dimensional,  or  planar  arrays,  are  central  to  the  focus  or  this  work,  we 
extend  further  the  random  linear  array  to  the  random  planar  array.  In  the  same  manner  as 
with  the  linear  arrays,  the  inter-element  spacing  is  no  longer  deterministic,  but  random  in 
both  coordinates.  In  the  previous  discussion  about  periodic  planar  arrays,  the  array  size 
contained  N  x  M  total  elements,  but  now  we  express  the  total  number  of  elements  as 
simply  N.  The  double  sum  in  Equation  (8)  is  merely  a  convenient  way  to  express  the 
sum  over  the  two  dimensions.  Here,  it  is  a  single  sum  because  we  explicitly  step  through 
the  coordinates  of  each  element  in  a  random  planar  array.  Thus,  for  N  total  elements,  the 
array  factor  becomes: 


AF(0,</>)=XV 


jk[(xn  sin  0  cos  (j)+ Px  )+(y/}  sin  0  sin  (j)+ (3  )] 


«=0 


(12) 


Substituting 


A,  =  -kxn  sinA,  cos  0O  and  =  -kyn  sin 6o  sin  0o 


(13) 


into  Equation  (12)  gives  the  array  factor  for  the  N  element  random  planar  array: 


AF(0,0)=XV 


jk[(xn  sinOcos(j)-xn  sin0o  cos Q0)+(yn  sin^sin^-y,,  sin0(,  sin0o)] 


(14) 


n—0 


where  (0O,0O)  determine  the  beam  pointing  angles.  This  equation  is  used  to  determine 
the  array  factor  in  all  subsequent  calculations  involving  random  planar  arrays. 


Random  arrays  exhibit  some  special  characteristics  that  differentiate  them 
from  periodic  arrays.  Some  details  of  random  arrays  are  presented  in  the  following 
subsections.  We  will  focus  our  attention  on  the  array  performance  characteristics  of 
grating  lobes,  the  average  array  factor,  main  lobe  and  beamwidth,  average  side  lobe  level 
and  peak  side  lobe  level. 


a.  Grating  Lobes 

Grating  lobes  are  a  result  of  the  periodic  spacing  of  the  elements  of  the 
array;  there  is  a  coherent  buildup  of  elemental  signals  at  angles  other  than  the  steering 
angle.  It  should  be  no  surprise  then  that  grating  lobes,  an  artifact  of  the  periodicity  of 
uniform  arrays,  are  significantly  reduced  in  random  arrays  because  of  the  random  nature 
of  the  elements  [19]. 
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b.  Main  Lobe  and  Beamwidth 

The  shape  of  the  main  lobe  (maximum  gain  and  HPBW)  and  the  shape  of 
nearby  side  lobes  are  not  significantly  affected  by  the  random  nature  of  the  elements  [19]. 
This  can  be  shown  by  examination  of  the  expression  for  the  average  power  pattern  of  a 
linear  random  array: 

(15) 

V  lM  J  JM 

where  |/0(#)|  is  the  design  array  power  pattern  and  N  is  the  number  of  elements.  More 
precisely,  near  the  main  beam  where  the  normalized  power  pattern  of  the  design  array  is 
approximately  one  (|/o(0)|"  ~1),  the  main  lobe  character  of  the  random  array  (maximum 

gain,  half  power  beam  width)  is  approximately  that  of  the  periodic  array.  Furthermore, 
the  variance  of  the  array  factor  is  given  by  [19]: 

<r»W  =  ^-(l-|/oW|2) 

where  f0(6 )  is  the  design  array  factor  and  N  is  the  number  of  elements.  As  the  number  of 
elements  N  grows,  in  the  limit  the  variance  approaches  zero: 

lim  al  =  0 

N—>°° 

This  simply  states  that  the  larger  the  number  of  elements,  the  smaller  the  variance,  and 
thus  the  closer  the  power  pattern  of  the  random  array  is  to  the  design  array  power  pattern. 

The  broadside  half  power  beam  width  (HPBW)  of  a  linear  geometric  array 
is  given  as  [18]: 

HPBW  =  0.886  — rad  =  50.8°  —  (16) 

L  L 

We  wish  to  compare  how  the  HPBW  of  a  random  array  compares  to  the  HPBW  of  a 
design  array  of  the  same  length. 

Figure  12  graphically  shows  how  Equations  (15)  and  (16)  compare  to  the 
power  pattern  and  HPBW  of  a  linear  array  -  the  design  array;  shown  here  for  arrays  of  5, 
8  and  20  elements  with  fixed  inter  element  spacing.  As  expected,  the  power  pattern  of 
the  random  array  is  almost  the  same  as  the  design  array  in  the  vicinity  of  the  main  beam, 
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but  forms  a  “pedestal”  of  magnitude  1/N  as  evidenced  in  Equation  (15).  The  random 
nature  of  the  array  will  cause  the  average  side  lobe  level  to  be  equal  to  the  pedestal.  We 
also  see  that  in  the  vicinity  of  the  main  lobe,  the  HPBW  is  approximately  equal  to  the 
design  array.  The  approximation  improves  as  the  number  of  elements  increases. 


- N=  8,  d!i=  0.5 

- Random  Army  Ave  PP 


0 

-5 

CD  '10 
■o 

-15 

^  -20 
-25 
-30 


0  (deg) 


Figure  12.  Comparison  of  the  Broadside  Power  Pattern  of  5,  8  and  20  Element  Linear 

Array  and  Average  Broadside  Power  Pattern  of  the  Same  Sized 
Random  Linear  Array 


The  results  of  Figure  12  are  important  because  they  give  us  an  envelope 
for  the  expected  power  level  of  the  random  arrays  to  be  developed  to  build  the  proposed 
communication  link  between  WSNs.  We  also  see  that  as  the  number  of  elements 
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increases,  the  HPBW  of  the  random  array  approaches  that  of  the  linear  array.  This  is 
important  as  we  can  then  easily  estimate  the  HPBW  of  a  random  array  based  on  using  a 
simple  formula  for  an  N  element  linear  array.  Therefore,  the  half  power  beam  width  of 
the  main  lobe,  for  either  a  random  array  or  a  uniform  linear  array  with  equal  spacing  and 
current  is  given  by  Equation  (16). 

c.  Average  Side  Lobe  Level 

While  the  main  lobe  occurs  in  the  direction  the  array  has  been  steered,  the 
side  lobes  are  those  that  occur  in  all  other  directions.  The  arithmetic  mean  of  the  peaks  of 
these  side  lobes  is  the  average  side  lobe  level.  In  contrast  to  the  structure  of  the  main 
lobe,  the  distant  side  lobe  structure  is  affected  by  randomness  of  the  element  spacing. 
This  can  be  seen  from  the  array  factor  in  Equation  (14),  where  xn  and  vn  are  random 
variables,  thus  the  phase  is  a  random  variable.  Off  of  the  main  beam,  where  the  power 

pattern  of  the  design  array  is  much  less  than  one  (|/o(0)|~  «1),  we  see  that  from  Equation 
(15)  the  power  pattern  of  the  random  array  is  dominated  by  the  1/N  term.  Since  the 
maximum  of  the  power  pattern  of  a  random  array  is  equal  to  that  of  the  linear  array,  and 
the  average  side  lobe  level  of  a  random  array  is  1  IN,  we  can  say  that  the  power  in  the 
average  side  lobe  level  is  equal  to  1/N  below  the  main  lobe  [19].  In  Figure  12  we  can  see 
this  graphically.  All  plots  are  normalized  to  0  dB.  In  Figure  12  we  see  the  average 
sidelobe  level  for  the  20  element  array  is  approximately  1/N  =  101og10(l  /  20)  =  -13  dB 
down  from  the  main  lobe. 

d.  Peak  Side  Lobe  Level 

While  the  maximum  of  all  side  lobes  of  a  periodic  array  is  generally 
termed  the  side  lobe  level  [18],  for  the  case  of  the  random  array  all  of  the  sidelobes  are 
random  in  nature  and  average  to  1/N.  The  peak  side  lobe  level  refers  to  the  maximum 
peak  sidelobe  value  relative  to  the  average  side  lobe  level  [19].  The  topic  is  not  easily 
summarized;  however,  a  rule  of  thumb  for  the  peak  side  lobe  level  is  that  it  is  unlikely  to 
exceed  the  average  side  lobe  level  by  10  dB  [19]. 


28 


In  conventional  periodic  arrays,  an  increase  in  performance  can  be 
achieved  by  using  a  non-uniform  current  distribution  across  the  elements  of  an  array.  For 
a  given  number  of  elements,  by  tapering  the  current  amplitude  near  the  edges  of  an  array, 
the  side  lobe  level  of  an  array  can  be  reduced  and  the  main  beam  width  increased  [18]. 
For  random  arrays,  however,  this  is  not  the  case;  a  tapered  current  distribution  has  only  a 
slight  impact  on  the  side  lobe  level  [19].  Therefore,  in  our  simulations  we  will  use 
uniform  current  amplitude  at  all  nodes  of  the  antenna  array. 

Having  discussed  the  characteristics  of  periodic  and  random  arrays,  we  now 
proceed  to  demonstrate  this  theory  using  a  MATLAB  simulation  of  the  random  planar 
array. 


B.  MATLAB  SIMULATION  RESULTS  -  RANDOM  PLANAR  ARRAY 

In  this  section,  the  periodic  planar  array  with  fixed  element  spacing  is  extended  to 
uniformly  distributed,  randomly  placed  elements  over  a  planar  area.  The  elements  are  z 
directed  half  wave  dipoles.  The  array  size  is  characterized  in  terms  of  wavelengths  and 
the  number  of  elements  N,  which  represents  the  total  number  of  elements.  This  is  in 
contrast  to  periodic  planar  arrays  which  commonly  list  the  number  of  elements  along 
each  axis,  e.g.,  Nx  =  5  and  Ny=  5  for  a  total  of  25  elements.  Here,  we  simply  use  N=  25 
to  indicate  that  there  are  25  elements  within  the  array. 

Elements  in  the  random  planar  arrays  are  restricted  to  being  no  closer  than  1/1 0th 
of  a  wavelength.  This  is  done  for  two  reasons,  to  better  model  a  physical  sensor  node  as 
it  is  unlikely  that  one  would  sit  on  top  of  another  and  also  to  reduce  the  influence  of 
mutual  coupling  between  antenna  elements.  Although  at  a  separation  distance  of  1/1 0th 
of  a  wavelength,  there  is  still  a  contribution  from  mutual  coupling  [18];  in  this  work  it  is 
ignored. 

In  Figure  14,  we  show  the  3  dimensional  gain  pattern  of  a  random  planar  array  of 
25  vertically  oriented  dipoles  with  the  elements  constrained  within  a  32.x  32,  area.  This 
plot  is  presented  to  demonstrate  the  pattern  of  the  random  planar  array  in  the  upper 
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hemisphere,  with  the  array  beam  steered  to  6=  90°  and  (j)  =  0°,  here  oriented  to  the  lower 
left.  Higher  gain  is  shown  in  bright  red,  while  lower  gain  is  shown  in  purple  and  blue.  In 
this  thesis,  we  are  interested  in  forming  the  terrestrially  bound  communication  link  and 
thus  we  will  focus  on  the  portion  of  the  pattern  that  lies  in  the  x-y  plane. 


Figure  13.  Three  Dimensional  Gain  Pattern  of  25  Element  3  Ax  3 A 
Random  Planar  Array  Steered  to  0=  90°  and  0=0°. 


In  the  following  sections,  through  simulation  we  examine  the  relationships 
between  the  physical  array  size  and  the  resulting  half  power  beamwidth;  the  gain  of  a 
fixed  size  array  and  the  number  of  elements  contained  within  it;  and  the  gain  of  the  array 
with  a  fixed  number  of  elements  versus  the  array  size.  These  results  will  be  useful  for 
developing  the  proposed  communications  link  between  WSNs  in  later  chapters.  Note 
that  since  the  WSNs  are  all  assumed  to  be  on  the  ground,  the  search  will  be  along  the 
horizon,  or  in  our  simulation’s  coordinate  system  (see  Figure  10)  we  have  6  =  90°,  and  <f b 
scanning  the  x-y  plane.  For  the  actual  sensors  (see  Figure  1)  the  transmit  and  receive 
antennas  are  vertically  oriented,  thus  we  only  simulate  the  vertical  polarization. 
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Additionally,  we  assume  here  that  the  radiation  efficiency,  77  is  one  so  that  the  terms 
directivity  and  gain  are  used  interchangeably. 

1.  Planar  Array  Size  versus  Half  Power  Beam  Width 

Given  a  random  planar  array  of  half-wave  dipoles  parallel  to  the  z-axis,  we  wish 
to  find  the  half  power  beam  width  in  the  direction  the  array  has  been  steered.  The 
broadside  half  power  beam  width  of  an  array  with  wavelength  land  length  L  is: 
1 

HPBW  =  50.8°  — .  This  also  applies  to  a  random  linear  array  of  length  L  [19].  We  now 
L 

consider  the  uniform  random  planar  array,  steered  to  0o  =  90°  and  0O  =  0°.  For 
observation  angles  in  the  x-y  plane,  0=  90°.  Here  we  are  interested  in  observations  angles 
near  the  main  beam,  i.e.  where  0  is  small  and  cos0  =  1 .  Thus  Equation  (14)  reduces  to: 

N- 1 

AF(0  =  90°, 0)  =  ^AnejkXn(cos*-l)ejkyn™* 

n= 0 

where  xn  and  vn  range  from  0  to  L  and  all  An  are  equivalent.  When  0  =  0°  this  is  then  the 
broadside  of  the  array  where  the  gain  is  maximum.  For  observations  where  0  is  small,  the 
length  L  in  the  y  dimension  determines  the  HPBW.  Because  cos0  changes  much  more 
slowly  than  sin0,  the  contribution  due  to  the  x  component  is  relatively  small.  Thus  the  x 
component  of  the  element  has  a  very  minor  influence  on  the  HPBW  in  the  broadside 
observation.  Therefore  the  expression  for  the  HPBW  of  a  linear  random  array  is 
approximately  valid  for  this  case  of  the  random  planar  array.  Because  the  distribution  of 
elements  within  the  Lx  L  area  is  random,  this  makes  the  choice  of  coordinate  system 
arbitrary.  If  the  pattern  is  constrained  to  the  x-y  plane,  the  expression  for  HPBW  is  valid 
for  an  arbitrary  steering  angle  in  the  x-y  plane.  For  larger  values  of  0,  the  array  factor  is 
dominated  by  the  random  character  of  the  elements  as  demonstrated  in  Figure  (12). 
Having  shown  that  the  HPBW  of  an  N  element  random  array  can  be  estimated  using 
Equation  (16),  we  now  validate  this  through  MATLAB  simulation  of  the  random  planar 
array.  In  this  simulation,  we  model  the  planar  array  by  arranging  elements  with  a 
uniformly  distributed  random  spacing.  In  this  simulation,  we  held  the  number  of 
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elements  constant  at  20.  We  then  varied  the  length  of  the  sides  of  the  array,  in  terms  of 
wavelength  and  calculated  the  HPBW  at  0  =  90°  and  0=0°  for  the  vertical  polarization. 
The  coordinate  system  is  the  same  as  in  Figure  10.  For  each  array  size,  we  generated  50 
realizations  of  the  array  of  random  elements  and  then  averaged  the  HPBW  values  from 
each  to  form  a  single  data  point.  At  each  array  size,  Equation  (16)  was  evaluated  and 
plotted  as  the  theory.  The  results  of  the  simulation  are  presented  in  Figure  14. 
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Figure  14. 


Broadside  Half  Power  Beam  Width  for  20  Element  Random  Planar  Array  and 
Theoretical  Half  Power  Beam  Width  Calculated  for  Increasing  Array  Size 


We  can  see  that  as  the  array  size  increases,  the  HPBW  becomes  smaller.  The 
results  show  good  agreement  with  the  theory,  indeed  the  analytical  expression  seems  to 
form  an  upper  bound.  From  this  result,  we  can  use  the  array  size  to  vary  the  beamwidth 
when  the  WSN  is  either  searching  for  other  WSNs  where  we  want  a  relatively  wide 
beam,  or  communicating  directly  with  another  WSN  where  we  want  a  high  gain  narrow 
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beam.  It  should  be  noted  that  as  the  array  size  grows  the  spacing  between  elements  will 
become  large  relative  to  a  wavelength,  especially  when  the  array  size  is  201  x  201 . 

2.  Planar  Array  Maximum  Gain  versus  Number  of  Elements  and 
Average  Side  Lobe  Level 

In  this  section  we  examine  how  the  maximum  gain  of  a  planar  array  is  affected  by 
increasing  the  number  of  elements  for  a  fixed  size  array.  In  Section  III. A. 2,  we  observed 
that  the  gain  of  an  A  element  random  linear  array  approximates  that  of  an  A  element 
linear  (design)  array;  additionally,  the  average  side  lobe  level  is  1/A  below  the  maximum 
gain  and  the  peak  side  lobe  level  is  no  more  than  10  dB  above  the  average  side  lobe  level. 
This  simulation  will  validate  these  observations  for  the  random  planar  array. 

In  Figure  15,  we  examine  the  maximum  gain  of  the  A  element  random  planar 
array  of  half  wave  dipoles  versus  the  number  of  elements.  The  physical  size  of  the  array 
is  held  constant  at  5  A  x  5  A  .  The  maximum  gain  is  calculated  for  the  steering  angle  0o  = 
90°  and  (j)0  =  0°,  and  is  plotted  in  blue.  Also  plotted  are  the  average  side  lobe  level  in  red 
and  peak  side  lobe  level  in  cyan.  The  average  side  lobe  level  was  computed  by  averaging 
all  the  peaks  of  the  side  lobes  for  each  given  number  of  elements.  For  each  number  of 
elements,  75  realizations  of  the  random  array  were  computed  in  order  to  form  an  average 
value  and  the  datum  plotted.  The  green  line  representing  the  theoretical  average  SLL  is 
from  Section  III.  A. 2. c  and  is  estimated  as  1/A  below  the  maximum  gain.  The  purple  line 
is  the  theoretical  gain  for  a  linear  array  of  A  dipoles  and  is  101og10(A)  plus  the  gain  for  a 
half  wavelength  dipole  of  2  dB  [18].  The  peak  side  lobe  level  in  cyan  is  presented  for 
reference.  As  discussed  in  Section  II.A.2.d  above,  the  peak  SLL  is  less  than  10  dB  above 
the  average  SLL.  The  cyan  line  shows  this  and  is  clearly  in  agreement  with  theory. 
More  exact,  but  significantly  more  complicated  predictions  for  the  peak  SLL  are 
presented  in  Chapter  8  of  [19]. 
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Figure  15.  Maximum  Gain,  Average  Side  Lobe  Level,  and  Peak  Side  Lobe  Level  as  a 
function  of  the  Number  of  Elements  in  a  Random  Planar  Array  of  Size  5  A  x  5A. 

Theoretical  Values  Are  Shown  for  the  Maximum  Gain  and  Average  SLL. 

In  Figure  15  we  see  that  the  maximum  gain  is  increasing  with  an  increasing 
number  of  elements  in  the  array  and  that  as  the  number  of  elements  increases,  the  average 
SLL  is  decreasing  farther  and  farther  below  the  maximum  gain.  The  peak  side  lobe  level 
is  also  decreasing  as  the  number  of  elements  increases.  This  is  in  agreement  with  theory. 

Note  that  the  maximum  gain  comes  from  the  main  lobe  of  the  array,  and  is 
therefore  unaffected  by  the  random  nature  of  the  array,  but  rather  is  proportional  to  N 
[18].  This  is  in  agreement  with  the  plot.  The  difference  in  gain  for  30  or  more  elements 
is  less  than  or  equal  to  1  dB  and  the  simulation  and  theoretical  results  converge  as  N 
increases.  Overall,  the  theoretical  and  simulation  curves  are  clearly  in  agreement. 
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From  the  theory  of  the  random  array,  the  average  side  lobe  level  should  be  less 
than  the  main  lobe  by  a  factor  of  1  IN.  This  is  depicted  by  the  red  and  green  lines  in  Figure 
15.  The  theoretical  and  simulation  results  again  show  excellent  agreement.  For  example, 
at  N  =  20,  the  average  SLL  should  be  13  dB  below  the  peak.  Here  the  peak  is 
approximately  17.0  dB  while  the  average  is  approximately  4.7,  yielding  a  difference  of 
12.3  dB.  At  N  =  50,  the  difference  is  16.3  dB  versus  a  calculated  17  dB. 

These  results  show  that  we  can  vary  the  amount  of  gain  for  a  fixed  size  array  by 
varying  the  number  of  elements  we  choose  to  include.  Thus,  for  a  particular  gain  level 
required  to  form  the  proposed  communications  link,  we  may  save  battery  power  by  only 
using  the  minimum  required  number  of  sensor  nodes  to  transmit  the  message  while 
simultaneously  extending  the  overall  lifetime  of  the  WSN. 

3.  Planar  Array:  Gain  versus  Array  Size  for  a  Fixed  Number  of 
Elements 

In  a  periodic  array,  gain  increases  with  array  size  because  elements  are  uniformly 
spaced.  In  a  WSN,  we  may  wish  to  achieve  a  particular  gain  level  without  increasing  the 
number  of  elements.  In  this  simulation,  we  demonstrate  that  we  can  thin  a  random  planar 
array  and  maintain  a  constant  gain.  We  fix  the  number  of  elements  in  the  array  at  20, 
then  increase  the  size  of  the  array  from  a  1A  x  1A  to  a  20A  x  20A  .  At  each  array  size,  we 
compute  the  maximum  gain,  average  SLL  and  peak  SLL  as  was  done  in  Figure  15. 

Figure  16  shows  the  results  of  this  simulation.  We  note  that  for  a  periodic  array, 
the  gain  is  actually  due  to  an  increase  in  the  number  of  elements.  Here  we  see  that  for  a 
fixed  number  of  elements  (20),  the  gain  shown  in  blue  is  relatively  constant  over  a 
varying  physical  size.  Since  the  number  of  elements  is  held  constant  and  the  physical  size 
is  increasing,  we  can  see  the  effect  of  reducing  the  element  density,  or  alternately, 
thinning  the  array.  Also  note  that  the  average  side  lobe  level  shown  in  red  is 
approximately  13  dB  below  the  main  lobe  gain,  which  is  given  by  Equation  (15)  and 
shown  in  green.  As  expected,  the  peak  SLL  does  not  exceed  the  average  SLL  by  more 
than  10  dB.  This  is  shown  in  cyan. 


35 


Figure  16.  Maximum  Gain,  Average  and  Peak  Side  Lobe  Levels  as  a  Function  of 
Random  Planar  Array  Size  for  Fixed  Number  of  Elements  =  20.  Theoretical 
Values  Are  Shown  for  the  Maximum  Gain  and  Average  SLL. 

As  in  Figure  15,  the  measured  maximum  gain  is  in  agreement  with  theory; 
however,  here  it  stabilizes  approximately  2  dB  above  the  theoretical  level  for  a  20 
element  array  of  dipoles  of  15  dB,  shown  in  purple.  The  simulation  appears  to  break 
down  for  array  sizes  less  than  4  wavelengths,  yet  it  quickly  converges  to  a  constant  gain. 

The  independence  of  gain  versus  array  size  is  an  important  result.  For  the  WSN, 
we  may  require  a  particular  gain  level  as  determined  by  the  number  of  elements  N 
composing  the  array.  From  the  results  of  Figure  16,  we  may  choose  any  N  elements  as 
long  as  their  distribution  is  still  random  within  the  bounded  area  [19]. 

4.  Summary  of  Results 

Figures  14,  15  and  16  demonstrate  some  important  results  that  we  can  employ 
when  building  the  proposed  communications  link  between  WSNs.  The  particular  WSNs 
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are  assumed  to  be  fixed  in  location,  but  the  inter-element  spacing  follows  a  uniform 
random  distribution  within  a  defined  boundary.  In  order  to  link  the  WSNs,  we  wish  to 
control  the  width,  gain  and  direction  of  the  array  beam.  We  have  shown  in  Equation  (14) 
that  we  can  steer  the  beam  in  an  arbitrary  direction.  In  Figures  14,  15  and  16  we  have 
demonstrated  control  over  the  beamwidth  and  gain. 

To  achieve  a  given  beamwidth,  we  choose  a  physical  size,  likely  a  sub-section  of 
the  WSN  under  consideration.  The  physical  size  controls  the  HPBW.  Figure  14  and 
Equation  (16)  demonstrate  this  behavior.  Within  this  physical  area,  we  have  control  over 
the  gain  as  demonstrated  by  Figure  15  by  choosing  the  number  of  elements  employed  in 
the  array. 

Alternately,  for  a  given  required  gain  level,  Figure  16  shows  that  we  can  thin  the 
array  (or  reduce  density)  dramatically  without  significant  loss  of  gain  or  the  effects  of 
grating  lobes.  It  does  not  matter  which  particular  sensor  nodes  within  the  WSN  are  used 
as  long  as,  on  average,  they  are  randomly  distributed  within  the  physical  area  chosen. 
This  is  important  as  we  may  wish  to  distribute  the  burden  of  transmitting,  and  thus 
consuming  energy  among  the  nodes  within  the  WSN. 

To  continue  from  Chapter  II,  in  the  first  example  using  a  pair  of  50  element 
random  planar  arrays,  with  transmitters  operating  at  full  power  (lmW  each)  and  a  margin 
of  15  dB,  we  determined  the  combined  gain  to  be  (G,  +  Gr)  =  15  -  Nt  +  44  dB.  From 
Figure  15,  we  see  that  by  choosing  A,  =  25,  we  have  a  transmitter  array  gain  Gt  of 
approximately  17.5  dB.  Therefore,  we  require  a  receiver  array  gain  G,  of  approximately 
16.5  dB.  This  can  be  achieved  with  17  receive  elements.  In  the  second  link  analysis 
where  the  transmitters  are  operating  a  minimum  power,  the  required  combined  gain  must 
satisfy  (G,  +  Gr)  =  15  -  N,  +  68  dB.  Again,  using  Figure  15,  we  see  that  using  45 
transmitter  elements  gives  us  a  transmitter  array  gain  Gt  of  19  dB.  The  required  receiver 
array  gain  G,  is  approximately  19  dB.  This  can  be  achieved  with  approximately  45 
receive  elements. 

Also  consider  a  20  A  x  20  A  area  with  240  elements  contained  within  the 
boundary.  Suppose  that  the  gain  requirement  is  such  that  we  only  need  15  elements.  We 
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are  still  free  to  choose  a  beamwidth  based  upon  the  communication  link  needs.  Suppose 
that  a  beamwidth  of  10°  is  required,  which  by  Equation  (16)  implies  a  5Ax5A  sized 
array.  Clearly,  the  WSN  may  subdivide  itself  into  5  A  x  5  A  segments  that  contain  at  least 
15  elements.  The  sub-sections  could  then  time  share  the  burden  of  transmitting.  Vincent 
in  [5],  provides  methods  for  power  conservation  and  burden  sharing  of  the  WSN. 

Armed  with  these  results  for  detennining  the  beamwidth  and  gain,  we  are  now  in 
a  position  to  apply  them  to  simulate  a  WSN  searching  for  other  WSNs. 


C.  SEARCH  METHOD 

In  order  for  one  WSN  to  find  another,  a  method  must  be  employed  to 
systematically  search  the  horizon.  In  this  thesis  we  have  chosen  to  employ  a  “lighthouse” 
approach  in  that  we  form  a  narrow  beam  and  transmit  in  a  given  direction,  then  steer  the 
beam  across  the  horizon  in  search  of  an  adjacent  WSN.  This  section  details  some  aspects 
of  this  method,  such  as  beam  steering,  allowing  for  response  time,  power  consumption, 
range  capability  of  the  array,  calculation  of  the  relative  range  and  bearing  between 
WSNs.  The  content  of  the  initiating  message  is  also  discussed  here. 

1.  Beam  Steering 

In  this  chapter,  we  have  established  that  a  random  planar  array  can  be  used  to 
form  a  beam,  and  this  beam  can  be  steered  to  a  given  direction.  Equations  (13)  and  (14) 
allow  us  to  steer  the  array  and  the  beamwidth  is  given  by  Equation  (16).  We  know  that 
beamwidth  is  related  to  the  size  of  the  chosen  array  (the  longer  the  baseline,  the  thinner 
our  beamwidth  and  vice  versa),  and  that  gain  is  related  to  the  number  of  elements  used. 
For  example,  to  scan  the  horizon  we  may  choose  a  half  power  beamwidth  of  30°,  transmit 
a  signal  in  search  of  another  WSN,  then  steer  the  beam  ahead  30°  and  transmit  again,  etc., 
until  we  have  moved  around  the  entire  horizon. 
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2. 


Look  -  Wait  -  Rotate  Beam 


This  search  method  is  such  that  the  searching  WSN  will  pick  a  start  direction, 
transmit  a  signal  in  that  direction  and  wait  for  a  response.  If  no  response  is  received,  then 
the  searching  WSN  will  advance  the  beam  and  “look”  in  the  next  direction. 

Since  each  node  in  the  WSN  acts  as  one  element  of  the  smart  antenna,  it  has  a 
particular  phase  or  weight  associated  with  it  to  steer  the  beam  in  each  look  direction.  The 
weights  can  be  stored  in  a  table  and  sequenced  through  in  order  to  advance  the  scan 
direction.  These  steering  weights  are  the  f3x  and  /3  of  Equation  (13).  Since  the 

geometry  of  the  WSN  is  fixed,  the  steering  weights  associated  with  each  sensor  node  (or 
transmit  element)  do  not  change  and  therefore  do  not  need  to  be  continuously 
recalculated.  All  nodes  within  the  WSN  transmit  with  the  same  amplitude.  As  noted  in 
Section  III.A.2.d  on  random  array  theory,  tapering  the  current  amplitude  weights  An  of 
Equation  (12)  with  respect  to  location  in  the  array  has  little  affect  on  the  performance  of 
the  array.  Since  all  elements  transmit  with  the  same  amplitude,  no  complex  current 
amplitude  weighting  scheme  need  be  designed,  simplifying  the  beam  steering. 

The  technique  of  rotating  the  beam  and  only  looking  in  one  direction  at  a  time 
was  chosen  to  avoid  having  several  WSNs  respond  at  once  as  might  occur  if  a  single 
omnidirectional  beacon  were  used.  If  a  single  beacon  signal  were  sent  omnidirectionally, 
the  searching  WSN  would  need  to  determine  the  direction  of  arrival  from  each  of  the 
responding  signals,  possibly  at  the  same  time. 

3.  Power  Conservation 

Power  conservation  is  not  a  major  focus  in  this  thesis;  however,  the  results  of  this 
work  are  meant  to  complement  other  research  on  energy  conservation  in  WSNs  [5],  [9]  , 
[21]. 

When  initially  searching  for  another  WSN,  the  array  does  not  know  the  distance 
to  the  other  array,  and  therefore  must  expect  the  longest  distance  possible.  This  implies 
using  all  elements  available  within  a  subset  of  a  given  area  to  form  the  initial  array,  to 
maximize  gain.  Since  we  initially  want  a  relatively  wide  beam,  for  fewer  total  beams  in 
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the  search  and  thus  fewer  transmissions,  these  constraints  limit  the  physical  size  of  the 
employed  sensor  transmitters.  Thus  a  tradeoff  will  be  required  to  optimize  power 
consumed  to  find  another  WSN  and  establish  the  link. 


4.  Range  and  Bearing 

It  is  assumed  that  one  goal  of  the  WSN  is  to  extend  its  lifetime  as  long  as  feasible. 
Since  the  sensor  nodes  are  battery  powered,  and  using  the  radio  transmitter  consumes  a 
relatively  large  amount  of  power  [5],  we  wish  to  minimize  the  number  of  transmitters 
used  in  forming  the  communications  link.  We  wish  to  only  use  the  minimum  number  of 
transmitters  required  for  a  given  range.  Since  a  higher  number  of  nodes  equates  to  a 
higher  transmission  gain  Gt,  from  the  Friis  equation,  we  can  determine  the  amount  of 
gain  required  for  a  given  distance  d  [12]: 

(  A  Y 

Pr=  —  P,G,Gr  (17) 

\Ajid  y 

where  P,  is  the  received  power,  G,  is  the  receiver  gain,  A  is  the  wavelength  and  d  is  the 
range.  From  Equation  (17),  we  see  that  a  higher  transmitter  gain  produces  a  higher 
received  power  Pr,  thus  for  a  given  receiver  gain  a  pair  of  WSNs  will  be  able  to 
communicate  at  a  longer  range  if  more  nodes  are  used. 


It  is  assumed  that  within  the  initial  exchange  of  infonnation  between  the 
searching  WSN  and  receiving  WSN  are  the  respective  latitude  and  longitudes.  From  this 
pair  of  lat/long  coordinates,  each  WSN  must  calculate  a  bearing  and  distance  to  the  other. 
For  a  given  pair  of  latitude  and  longitude  coordinates  ( lat1,longl )  and  (lat2,long2) ,  the 

distance  between  them  is  given  by  the  Haversine  formulas  [22],  Here,  the  WSNs  are 
assumed  to  be  on  the  order  of  a  kilometer  apart  and  the  flat  Earth  approximation  to  the 
Haversine  formulas  is  applicable  [22].  The  separation  distance  D  (in  km)  and  bearing 
<D  (in  radians)  are  calculated  using  Equations  (18)  and  (19): 


O  =  tan  1 


^  sin(A/ong)cos(/at2) 

cos(toj )  sin(/a?2 )  -  sin(to, )  cos  {lat2 )  cos  ( A  long ) 


(18) 

(19) 
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where  a  =  —  —  latl ,  b  =  —  -  lat2  and  A  long  =  long2  -  Iongl .  A  simple  java-script  of  the 

Haversine  formulas  for  the  flat  Earth  assumption  is  available  at  [23].  Once  the  range  and 
bearing  are  known,  the  WSNs  can  calculate  the  required  gain  to  establish  the  link.  This 
gain  calculation  was  the  subject  of  Chapter  II.  For  the  sensor  nodes,  presumably 
deployed  on  the  ground,  the  antenna  height  hant  is  approximately  0.1  m.  The  radio 
horizon  for  this  node  is  cl  —  4A2yJhant  [24]  or  1.3  km.  Thus  the  flat  Earth  approximation 
used  in  Equations  (18)  and  (19)  for  this  short  distance  seems  applicable. 

We  now  demonstrate  a  specific  simulation  of  a  WSN  forming  a  beam  and 
scanning  it  across  the  horizon. 


D.  MATLAB  SIMULATION  RESULTS  -  SCANNING  BEAM 

We  now  wish  to  combine  the  concepts  of  this  chapter  on  random  planar  arrays 
with  the  idea  of  a  searching  WSN  that  transmits  a  beam  in  a  rotating  fashion  across  the 
horizon.  We  begin  by  choosing  a  half  power  beamwidth  of  approximately  20  degrees, 
from  Equation  (16),  and  set  the  array  size  at  2.5  A  x  2.5  A  .  In  this  example  we  do  not 
have  a  specific  gain  requirement  to  meet,  so  we  pick  40  elements  to  obtain  a  gain  of  16 
dB  as  a  convenient  number.  To  steer  the  beam,  we  change  the  phase  of  each  element  in 
accordance  with  Equation  (13). 

Figure  17a  shows  the  realization  of  a  random  array  of  40  half  wave  dipole 
elements  distributed  over  a  2.5Ax  2.5A  area,  and  in  Figure  17b  the  coordinate  system  is 
presented  for  reference.  The  minimum  spacing  between  elements  is  set  to  0.  IX,  but  would 
need  to  be  increased  to  realize  this  example  with  the  MICAz  mote.  The  0.1X  spacing 
found  in  some  parts  of  this  realization  is  much  denser  than  a  traditional  planar  array  of 
the  same  area.  As  before,  is  used  to  visualize  the  coordinate  system  by 

representing  an  arbitrary  angle.  The  cyan  circle  depicts  the  locus  of  observation  points, 
beginning  at  9  =  90°  and  with  (j)  =  0°  on  the  x-axis.  Figure  18  shows  the  gain  pattern  of 
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the  array  steered  to  0o= 90°,  (f>0  =  0°.  Figure  19  shows  polar  plots  of  the  gain  patterns  for 
the  first  six  beams  to  envision  the  scanning  perfonned  by  the  array. 
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Figure  17.  Searching  WSN: 

(a)  Realization  of  a  40  Element,  2. 5 A  x  2.5/1  Random  Planar  Array 

(b)  Coordinate  System  with  Scan  Plane  at  9=  90°  Shown  in  Cyan 


In  Figure  18,  we  see  the  main  lobe  with  a  peak  gain  of  16.7  dB,  comparable  to  the 
theoretical  gain  of  16  dB  for  a  40  element  array.  A  random  shape  characterizes  the  side 
lobes,  with  the  expected  pedestal  nature  shown  in  Figure  12.  The  average  side  lobe  level 
for  all  scan  angles  is  at  17.1  dB  below  the  main  lobe.  The  half  power  beam  width  in 
Figure  18  is  22°,  in  close  agreement  with  the  expected  HPBW  of  50.8°/2.5  =  20°. 
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Figure  18.  H-plane  Gain  Pattern  for  a  40  Element  Random  Planar  Diploe  Array 

with  Beam  Steering  Angle  0o=9O°,  (j)0  =  0° 


Figures  19a-f  show  the  gain  pattern  for  six  consecutive  beams  in  a  polar  plot 
format.  Figure  19a  is  the  same  data  as  Figure  18  repeated  in  a  polar  format. 
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Gain  Pattern,  Scan  Angle  =  0.0 


Gain  Pattern,  Scan  Angle  =  80.0 
105-  90‘  75- 


Gain  Pattern,  Scan  Angle  =  40.0 


Gain  Pattern,  Scan  Angle  =  1 20.0 
105-  90‘  75- 


Gain  Pattern,  Scan  Angle  =  1 60.0 
105-  90-  75- 


Gain  Pattern,  Scan  Angle  =  200.0 
105-  90-  75- 


Figure  19.  H-plane  Patterns  for  a  Random  Planar  Array  of  40  Dipoles  as  Beam  is 

Scanned  From  0  to  200°  in  40°  Steps 
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We  see  that  as  the  beam  is  steered  in  each  direction,  it  retains  approximately  the 
same  gain  and  average  side  lobe  levels  on  each  plot. 

It  should  be  noted  that  in  this  example  we  used  a  20°  beamwidth  to  reduce  the 
number  of  beams  required  to  scan  the  entire  horizon,  and  thus  the  total  power  to  complete 
one  sweep.  However,  to  obtain  this  beamwidth,  the  array  size  is  limited  to  2.5/1  x  2. 5 A, 
which  in  turn  limits  the  total  number  of  elements  that  can  physically  fit  within  that  area. 
A  reduced  number  of  elements  results  in  a  lower  peak  gain.  Thus  a  tradeoff  is  required 
between  the  total  number  of  elements  (peak  gain)  and  the  beamwidth  (size  of  the  array 
and  total  energy  consumed  to  scan  the  horizon). 

These  plots  demonstrate  that  the  beam  of  a  WSN  composed  of  uniformly 
distributed,  randomly  placed  sensor  nodes  can  be  used  to  form  a  beam  and  rotate  that 
beam  across  the  horizon.  The  characteristics  of  the  antenna  beam  are  consistent  with  the 
theory  presented  in  this  chapter  on  random  antenna  arrays.  This  completes  the  discussion 
on  the  transmitting  or  searching  WSN.  As  the  beam  from  the  transmitting  array  scans  the 
horizon,  it  is  assumed  to  illuminate  a  receiving  WSN,  which  must  recognize  and  reply  to 
this  signal. 

In  this  chapter,  we  presented  results  demonstrating  control  over  the  half  power 
beam  width,  and  gain  and  direction  of  the  radiation  pattern  for  a  random  planar  array. 
We  demonstrated  that  we  could  select  a  wide  beam  for  the  search  pattern  and 
simultaneously  select  the  gain.  This  beam  was  steered  in  a  lighthouse  fashion  simulating 
the  search  for  and  hopefully  illumination  of  an  adjacent  WSN. 

We  now  turn  our  attention  to  the  receiving  array. 
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IV.  DIRECTION  OF  ARRIVAL 


The  approach  of  this  chapter  is  essentially  from  the  viewpoint  of  the  receiving 
array.  When  the  receiving  WSN  detects  a  signal  from  the  searching  WSN,  it  needs  to 
know  the  relative  geometries  of  the  two  WSNS  in  order  to  form  its  array  and  then  send  a 
reply.  The  location  information  of  the  sending  WSN  is  contained  within  the  message 
signal  transmitted,  but  the  direction  of  arrival  (DOA)  must  be  calculated.  This  chapter 
establishes  the  need  to  blindly  determine  the  DOA  of  the  received  signal  and  presents  a 
short  survey  of  techniques  available  in  the  literature.  A  spread  spectrum  space  division 
multiple  access  technique  is  introduced  and  demonstrated  first  with  a  linear  array,  and 
then  applied  to  the  random  planar  array.  A  simulation  of  the  message  extraction  from  a 
received  signal  is  presented.  The  strengths  of  this  technique  as  well  as  the  importance  of 
the  spread  spectrum  pseudorandom  noise  codes  are  explained. 

A.  DETERMINATION  OF  THE  RELATIVE  GEOMETRIES 

Here  we  consider  two  of  the  WSNs  depicted  in  Figure  3,  WSNs  “A”  and  “B.”  We 
consider  that  WSN  A  has  sensed  a  phenomenon  and  needs  to  pass  this  data  to  the  next 
WSN,  say  WSN  B.  In  Figure  20,  WSN  B  has  been  illuminated  by  the  searching  WSN  A 
and  now  needs  to  determine  the  DOA  of  the  signal.  The  receiving  WSN  B  will  form  a 
random  planar  array  using  its  sensor  nodes  to  analyze  the  incoming  signal  and  determine 


the  DOA. 


WSN  A 


WSN  B 


Figure  20.  Signal  From  WSN  A  Illuminating  WSN  B 
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The  receiving  WSN  will  then  use  this  calculated  DOA  and  the  location 
information  encoded  within  the  transmitted  signal  to  determine  the  relative  geometries 
and  form  the  optimal  return  beam  to  the  searching  WSN. 

1.  Link  Initiating  Signal 

In  Figure  20,  the  receiving  array  has  no  prior  knowledge  of  the  searching  WSN 
from  which  the  initiating  signal  is  coming.  It  will  have  to  determine  the  direction  of 
arrival  by  processing  the  signals  received  at  each  of  the  sensor  nodes.  Section  B  of  this 
chapter  discusses  this  in  detail. 

Some  sensor  nodes  have  the  capability  to  operate  in  a  power  saving  or  “sleeping” 
mode.  This  thesis  does  not  address  the  sleeping  or  transition  of  modes  for  each  of  the 
sensor  nodes.  Instead,  we  assume  that  either  the  search  beacon  contains  a  “wake  up” 
signal  or  some  other  mechanism  has  allowed  all  sensors  in  the  receiving  array  to  be 
turned  on  for  receiving  signals. 

2.  Return  Signal 

We  assume  that  in  addition  to  the  location  of  the  sending  array,  the  message  in  the 
search  signal  contains  the  “name”  or  “address”  of  the  searching  WSN,  e.g.,  its  MAC 
address  and  its  lat/long  coordinates.  The  receiving  WSN  is  then  in  a  position  to  solve  its 
half  of  the  link  equation  and  form  a  return  signal.  For  example  it  may  have  determined 
that  WSN  A  at  bearing  040°  and  range  500  m  has  requested  a  link.  When  it  forms  its 
reply  or  return  message  to  the  originator,  it  must  include  its  lat/long  coordinates  as  part  of 
the  initial  handshake.  The  first  task,  however,  is  to  determine  the  DOA  of  the  query. 


B.  DOA  DETERMINATION 

Godara  [10]  provides  a  summary  and  comparison  of  various  direction  of  arrival 
methods  including  spectral  estimation,  minimum  variance  distortionless  response 
estimator,  maximum  likelihood,  maximum  entropy,  least  mean  squares,  conjugate 
gradient,  Multiple  Signal  Classification  (MUSIC)  and  Estimation  of  Signal  Parameters 
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via  Rotational  Invariance  Technique  (ESPRIT).  MUSIC  and  ESPRIT  are  also  known  as 
eigenstructure  or  subspace  methods  and  can  be  computationally  intensive  [10].  ESPRIT 
and  most  versions  of  MUSIC  require  periodic  arrays  [20],  thus  making  them  unsuitable 
for  the  WSN  in  our  problem.  Time  difference  of  arrival  (TDOA)  is  another  technique 
used  in  many  applications  [9].  The  direction  of  arrival  problem  is  also  known  as  angle  of 
arrival,  spectral  estimation,  or  bearing  estimation. 

1,  The  Requirement  for  a  Blind  DOA  Algorithm 

From  the  perspective  of  the  receiving  WSN,  when  a  signal  is  received  from  the 
searching  WSN,  its  DOA  is  unknown.  The  receiving  WSN  must,  therefore,  determine 
the  direction  from  which  the  searching  WSN  is  transmitting  without  the  aid  of  a  reference 
signal  as  is  used  by  some  DOA  methods.  DOA  determination  methods  that  do  not  use 
prior  knowledge  of  signal  direction  or  some  sort  of  reference  signal  to  form  an  initial 
estimate  for  the  DOA  are  tenned  blind  DOA  methods  [10]. 

A  new  method  based  upon  spread  spectrum  techniques  for  space  division  multiple 
access  (SDMA)  applications  that  does  not  involve  iterative  matrix  solutions,  reference 
beams  or  specific  array  geometry  has  been  presented  and  described  in  a  patent  by  Elam 
[25].  Because  this  method  does  not  possess  the  computational  intensity  of  the  above 
methods,  is  suitable  for  random  arrays  and  has  the  capability  for  blind  DOA 
determination,  we  chose  this  method  for  the  receiving  array.  The  rest  of  this  chapter 
discusses  and  applies  the  SDMA  technique  to  our  problem. 


C.  SPACE  DIVISION  MULTIPLE  ACCESS  RECEIVER 

The  direct  sequence  spread  spectrum  (DSSS)  technique  proposed  by  Elam  [25] 
and  as  described  by  Gross  [20]  is  summarized  in  this  section.  Analogous  to  the  code 
division  multiple  access  (CDMA)  receiver,  this  space  division  multiple  access  receiver 
(SDMA)  spatially  separates  received  signals  using  DSSS  techniques.  The  method 
involves  chipping  the  phase  of  the  received  signals  at  each  array  element  with  individual 
spreading  sequences.  In  computer  memory  identically  chipped  virtual  signals  are 
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generated  from  a  virtual  receiver  array  for  /  expected  directions  of  arrival.  The  received 
signals  are  then  summed  and  a  quadrature  correlation  is  performed  between  the  received 
signal  and  each  of  the  I  expected  DOA  virtual  signals.  Correlation  values  (R,)  exceeding 
a  threshold  are  identified  as  signals  received  from  an  expected  DOA;.  From  the  phase  of 
the  correlation  Rt,  the  message  signal  information  may  be  extracted.  Interfering  signals 
are  not  well  correlated  and  thus  minimized  [20].  This  is  similar  in  concept  to  a  matched 
filter  bank  in  a  code  division  multiplexing  receiver. 

The  time  sample  of  the  received  signals  is  of  length  T.  With  this  time  sample,  M 
pseudorandom  chips  are  applied  to  the  received  signal  at  each  receive  element.  Each 
receive  element  has  a  separate  pseudorandom  noise  (pn)  spreading  code.  The  spreading 
sequences  can  be  thought  of  as  pseudorandom  array  weights.  If  viewed  in  the  context  of 
antenna  array  weights,  during  the  time  period  T,  the  array  has  been  “steered”  by  M  sets  of 
pseudo-random  weights.  Over  all  M  chips,  the  average  pattern  approaches  that  of  an 
omni-directional  array  [20].  Viewed  another  way,  the  “look  direction”  beams  are  created 
through  the  correlation  with  the  virtual  beams.  Thus,  steering  the  “look”  direction  of  the 
physical  array  is  done  through  the  software  that  creates  the  virtual  array  expected  DOAs. 


We  now  discuss  the  receiver  in  greater  detail.  Figure  21  illustrates  the  receiver 
(after  [20,25]).  In  the  following  sections  we  examine  the  three  main  components  of  the 
receiver:  the  actual  array,  the  virtual  array  and  the  correlation  and  phase  detection.  We 
then  describe  the  method  for  extracting  the  message  in  the  received  signal. 
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SDMA  Receiver  (after  Gross  and  Elam) 
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Figure  21. 


1. 


The  Actual  Array 


The  actual  array  represents  the  physical  receive  antenna  elements  in  the  receiving 
WSN  with  a  given  geometry  and  L  signals  arriving  at  the  array.  The  signals  are  assumed 
to  have  been  phase  shift  keyed  (PSK).  The  L  phase  modulated  incident  signals  s,(t) 

arrive  at  the  receiver  and  form  the  received  signal  xn(t)  at  each  of  the  A  antenna  elements. 
The  A  elements  need  not  be  uniformly  spaced,  but  can  be  arbitrarily  located.  Indeed, 
randomly  placed  elements,  in  a  two  or  three-dimensional  array  helps  the  performance  by 
producing  independent  phasing  at  each  element  for  each  expected  angle  of  arrival.  Each 
received  signal  xn(t)  is  then  spread  or  chipped  with  a  pseudorandom  noise  (pn)  code 


associated  with  that  antenna  element,  producing  the  spread  signal  rn  (?) ,  and  finally  the 
summed  signal  r(t) .  The  signal  is  sampled  over  a  time  period  of  t  to  t+T.  The  received 
signal  at  each  array  element  can  be  expressed  as: 

L 

Xn(t)='Za„,lS,(t)  (20) 

1= 1 

where  for  a  uniform  linear  array  with  k  as  the  wavenumber  and  d  the  distance  between 
elements,  the  array  factor  contribution  an ,  for  element  n  from  signal  /  is  given  by: 


n  —  jkd(n-l)  sin  dt 
Un,l  ~  e 

and  the  baseband  phase  modulated  received  signal  of  interest  s,  (/)  is  given  by: 


(21) 


s,(t)=  ejQ(,)  (22) 

From  Equation  (20),  we  can  express  the  received  signal  at  each  of  the  receive  elements  as 
an  A  x  1  column  vector: 


x(t)  =  [xl(t)  x2(t )  ...  xN(t)J 

For  A  elements  and  L  signals,  the  received  signal  can  be  expressed  in  matrix  form  as: 

x(t)  =  As(t)  (23) 

where  A  is  the  NxL  matrix  of  the  steering  vector  elements,  and  .?(/)  is  an  Ax  1 
column  vector  of  the  incident  signals  S/(t) .  The  matrix  A  can  be  compactly  expressed  as 


a  1  x  L  row  vector  of  the  steering  vectors  a[  (0, ) : 
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x(t)  =  \a[(G i)  d'2(02 )  •••  a[(0j] 


(24) 


e?2(0 

/W) 


or  A  can  be  expressed  showing  the  expanded  steering  vectors  in  order  to  see  the 
individual  contribution  of  the  nth  element  and  /th  signal: 


l 

1 

1 

1 

1 _ 

II 

/^s 

gjkd  sin  0, 

^jkd  sin  02 

gjkd  sin0L 

^(JV-Osinfl, 

gjkd(N-V)s\n0L 

i 

1 _ 

For  example,  for  L=  1,  x(t )  is  given  as: 


x(t)  = 


gtfitOgjfaisine, 


gtfi  (t)gjkd(N-\  )sine, 

The  pn  chipped  signal  r(t)  is  described  in  matrix  form  as: 


(25) 


(26) 


where 


r(t)  =w(t)x(t) 


(27) 


w(t)  =  [w0(0  w,  (0  •••  wA,_,  (f)] 


(28) 


is  the  1  x  N  row  vector  of  the  spreading  codes.  These  spreading  codes  wn(t)  will  be  used 
to  phase  modulate  the  array  outputs  xn(t).  Each  vector  wn(t)  contains  M  chips  of  length  x, 
chosen  such  that  the  length  of  the  vector  is  Mr  —  T ,  the  sample  time.  The  number  of 
chips  M  should  be  chosen  such  that  it  is  greater  than  the  Nyquist  rate  of  the  received 
baseband  signal  modulation  [20].  This  implies  the  phase  modulation  of  the  received 
baseband  signal  is  nearly  constant  over  T.  The  time-bandwidth  product  of  the  sample 
time  T  and  the  message  signal  bandwidth  B  should  be  TB  <  0.25  [20]. 


In  spread  spectrum  applications,  there  are  many  types  of  pseudorandom  noise 
spreading  codes  available  including  maximal-length  sequences  or  ///-sequences,  Gold 
codes  [26],  and  Golay  codes  [27].  Walsh-Hadamard  codes  [27]  and  Welti  codes  [28]  can 
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be  used  as  orthogonal  covering  spreading  codes.  In  our  implementation  we  have  chosen 
Walsh-Hadamard  codes  because  of  the  simplicity  of  implementation. 

The  cross  correlation  of  codes  within  the  same  family  provides  one  measure  of 
how  well  a  spreading  code  works  [12].  Codes  that  have  a  minimum  cross-correlation 
will  cause  a  minimum  interference  with  other  codes  because  signals  spread  with  one  code 
will  not  correlate  with  another  code  and  be  suppressed  [26].  The  orthogonality  of  the 
Walsh-Hadamard  or  Welti  codes  provides  separation  for  each  of  the  channels,  here 
represented  by  receive  elements  which  separate  out  the  signals  from  each  direction  [20]. 
Welti  codes  have  been  reported  to  produce  the  best  results  in  this  application  [29], 

Figure  22  shows  the  effect  of  the  spreading  codes  on  the  received  signal  at  one  of 
the  antenna  elements.  Here  we  show  the  received  signal  phase  £n(t)  at  one  of  the  antenna 
elements  in  blue,  the  Walsh-Hadamard  spreading  code  wn  in  red,  and  the  real  part  of  the 
chipped  signal  rn(t)  in  purple.  The  chipping  process  is  repeated  at  each  of  the  N  receive 
elements,  each  element  having  its  own  individual  Walsh-Hadamard  code. 
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time(t)  x  10^ 


Figure  22.  Received  Signal  Phase  C,( t),  Walsh-Hadamard  Spreading  Code  wn, 

and  the  Real  Part  of  the  Chipped  Signal  rn  (l  ) 


Key  to  this  technique  is  the  formation  of  the  virtual  array  of  signals  that  “look”  in 
the  expected  directions  of  arrival.  The  generation  of  the  virtual  signals  is  similar  to  that  of 
the  actual  array.  We  continue  with  the  description  of  the  SDMA  receiver  by  examining 
the  virtual  array  component. 

2.  The  Virtual  Array 

Figure  23  graphically  shows  this  component  of  the  receiver.  Within  computer 
memory,  a  “virtual  array”  is  configured  according  to  the  geometry  of  the  physical  array, 
containing  the  same  number  of  elements.  The  virtual  array  is  steered  to  an  expected  DOA 
by  computing  the  array  factor  for  that  expected  direction  6,.  The  expected  direction  Oi 
can  also  be  thought  of  as  the  beam  steering  angle  0o  .  The  array  factor  for  each  of  the  N 
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elements  is  chipped  with  the  identical  spreading  codes  as  the  physical  array.  This  is 
repeated  and  stored  for  each  of  the  /  expected  directions  of  arrival  ft.  The  /  expected 
directions  can  be  chosen  to  avoid  known  interfering  directions,  prior  knowledge  of  the 
incoming  signal,  multipath  [20,25],  etc.  The  expected  directions  cannot  be  chosen  any 
finer  than  physically  allowed  by  the  actual  array. 


Virtual  Array 


Figure  23.  Expected  Signal  Generation  by  the  Virtual  Array  in  Memory 


For  each  of  the  /  expected  directions,  we  can  express  the  summed  virtual  signal 
as: 


vi(t)=m%(Oi)  (29) 

where  the  steering  vector  for  each  expected  direction  is  given  by  the  N  x  1  column 
vector: 


1 

jkd  sin  ft 

eJ 


gjkd(N~  l)sin0, 

and  the  spreading  codes  are  the  same  as  before: 

w(t)  =  \w0{t)  Wjft)  •••  hv_i(0] 


(30) 


Having  formed  both  the  received  signal  and  the  array  of  virtual  signals,  the 
receiver  now  performs  a  correlation  computation  to  establish  the  DOA  of  signals  incident 
to  the  receiver. 
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3.  Signal  and  Phase  Correlator 

The  signal  and  phase  correlator  correlates  the  actual  received  signal  r(t)  with 
each  of  the  /  virtual  signals  v(.  it) ,  over  the  sample  time  T.  Figure  24  illustrates  the 
operation. 


lit) 


Figure  24.  Signal  and  Phase  Correlator 

The  correlation  computation  is  expressed  as: 

ft,  =  j;  r(t)  v*(t)  dt  (31) 

where  R,  is  the  complex  correlation  value  for  direction  0i .  The  magnitude  of  /?,  will  be 
plotted  for  each  angle  0,  becoming  the  array  factor  for  the  signals  present.  The  threshold 
detector  depicted  compares  the  magnitude  of  R,  against  a  predetermined  value  and 
declares  a  signal  present  or  not  at  this  angle.  If  a  signal  is  present,  the  data  is  extracted  as 
explained  in  the  next  subsection.  This  is  similar  to  a  matched  filter  or  correlation  detector 
operation  [27]. 

The  correlation  values,  when  the  pn  codes  for  the  actual  and  virtual  signals  match, 
become  the  array  factor  for  the  antenna  array  in  the  direction  of  the  virtual  signal.  To  see 
how  this  works,  we  examine  the  correlation  in  more  detail.  From  Equations  (22)  and  (27), 
the  actual  signal  received  at  the  antenna  array  r(t)  with  N  array  elements  is  expressed  as: 
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(32) 


r(t)  =  X 


ej(n~l)kd  sin@0  gXO)f  2  W"  ^ 


n= 1 


where  the  spreading  codes  wn  are  modulating  the  phase  of  the  received  signal.  The  signal 
generated  in  memory  for  the  i,h  expected  direction  and  modulated  with  the  identical 
spreading  codes  is  from  Equation  (29): 


N  ,n_ 

( t )  =  X me‘ e  lWm(  ’ 


(33) 


Substituting  Equations  (32)  and  (33)  into  (31)  and  noting  that  the  signal  has  been  sampled 
in  time  through  the  spreading  process  gives  the  expansion: 


R 


:=JX 


.n  JV 

ej(n-l)kdsmeo  eX(l)eJ2Wn(°^ 

m= 1 


e-Km-l)kdMi6i 


(34) 


Rearranging  the  element  summation  yields 


N  N 


*<=JXX 


j(n-Y)kd sin 6„  gXWf  2  M"  ^ e~j(jn-l)kd sing,  7  2  (0 


(35) 


^  «=1  m=l 


Since  the  array  factor  terms  are  independent  of  time,  they  can  be  moved  outside 
the  integral.  Recall  that  we  chose  the  time  sample  period  T  such  that  the  phase  £(t)  of  the 
signal  s(t)  was  approximately  constant  over  the  sample  period  and  can  be  approximated 
as  its  average  £  ;  thus,  it  can  be  moved  outside  the  integral  as  well: 


ft  _  ’y'  ’y'  ej(n-l)kdsmeoe-j(m-\)kdsmei  ^X  J  J 2W"<'^ g  7 


8(n-m)=- 


1  if  n=m 
0  if  n^m 


(36) 


The  integral  yields  a  delta  function  8{n-m )  because  wn  and  wm  are  orthogonal  [12]. 
Thus  Ri  reduces  to: 


ft  —  yy(»-l)M(rin0.-sin8,)  gX 


(37) 


Equation  (37)  is  equivalent  to  the  array  factor,  here  multiplied  by  the  constant 
phase  of  the  received  signal.  Consequently,  the  correlation  can  be  expressed  in  terms  of 
magnitude  and  phase  as: 

R,.=|RiyPi  (38) 
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where  /?,.  is  the  magnitude  of  the  complex  valued  correlation  and  p,  expresses  the  phase 

of  the  received  signal.  The  magnitude  of  the  correlation  contains  the  DOA  information 
and  the  phase  term  contains  the  message  information. 

4.  Message  Extraction 

The  phase  of  the  signal  s(t)  can  also  be  estimated,  and  the  message  data  extracted 
using  this  SDMA  method.  Over  the  time  sample  period  T,  the  phase  of  signal  s(t)  was 
assumed  to  be  approximately  constant.  Thus,  when  a  signal  is  detected  from  direction  i, 
the  phase  and  thus  the  message  can  be  approximated  as  p, .  The  correlation  integral 
averages  the  received  signal  phase  over  T  [20],  therefore 

—  1  rt+T 

A“  C  =  A  (,m  (39) 

T  Jt 

where  £  is  the  average  signal  modulation  during  the  time  period  T  at  ft.  By  calculating 
p,  at  each  angle  where  a  signal  has  been  determined  present,  we  can  simultaneously 
extract  more  than  one  message  signal  from  the  correlation  data.  The  process  is  described 
in  further  detail  in  [20,  25]. 

5.  The  Importance  of  the  Spreading  Codes 

One  may  ask  what  advantage  do  the  spreading  codes  provide  to  the  overall 
receiver?  The  answer,  in  short,  is  that  the  weight  functions  or  the  chipping  sequences 
provide  orthogonality  of  each  received  signal  to  all  of  the  incorrect  DOA  virtual  signals. 
To  examine  this,  we  consider  an  example  with  N  =  4  elements;  the  signal  is  incident  from 
direction  0inc  and  wn,  k  and  d  are  defined  as  before.  The  received  signal  is  then: 

r(t)  =  eKl(,)  [ w0(t)es°  +  w^e**™0*  +  Wi{t)eik2d^  +  w3(t)eJk3dsme~  ] 

Recalling  that  the  time  sample  period  T  was  chosen  such  that  s(t)  was 
approximately  constant  in  phase,  the  leading  exponential  can  be  suppressed  as  a  constant. 
If  r(t)  is  not  chipped,  the  time  sample  is  reduced  to  a  single  sample  point  in  time,  s(t0), 
and  the  expression  for  r(t)  becomes  the  summation: 

N  N 

KO  =  X<  =  S^("'1,Wsme"'c  (4°) 

n= 1  n= 1 
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where  dn  represents  an  array  elements’  contribution  to  the  array  factor. 

For  each  of  the  /  expected  directions,  the  virtual  signal  is  expressed  as: 

V,«„)=  (41) 

n= 1  n= 1 

where  d„i  represents  the  virtual  array  elements’  contribution  to  the  array  factor. 

Using  Equations  (40)  and  (41)  in  the  correlation  computation  yields: 

r, =/•«>;<(„)= £ix<c  =  'Ll <42) 

n= 1  m= 1  n= 1  m= 1 

This  expression  is  the  correlation  sum  that  results  when  the  spreading  codes  are  not  used. 
Now  substituting: 

ys inc  =  kd  sin 0 inc  and  i//(  =  kd  sin  (9,  (43) 

into  Equation  (45)  and  expanding  the  sums,  the  correlation  sum  takes  the  form: 

R  =  [eiOV,„  +  ejVmc  +  eJ2Vix  +  eJ3Vinc  J [e-J0¥l  +  e-M  +  e-j2Vi  +  e-j3Vi  j  (44) 

The  correlation  magnitude  should  be  at  a  maximum  when  the  virtual  angle  6 ,  matches 
the  incident  signal  angle  dinc,  i.e.,  y/!m:  =  i//(  =  kd sm  6 !nc=  a  constant.  After  some 
algebraic  manipulation,  Equation  (44)  can  be  written  as: 

3  1 

Rmm  =  4  -l-  —  cos  (kd  sin  Qinc )  +  cos(2  kd  sin  6inc )  +  -  cos(3  kd  sin  Qmc )  (45) 

While  our  intent  was  that  the  correlation  reached  its  maximum  when  the  expected 
angle  matched  the  actual  incident  angle,  this  clearly  is  not  the  case.  In  Equation  (45),  the 
cosine  and  sine  terms  will  influence  the  value  of  R,  depending  upon  the  particular  value 
of  0  for  the  incident  signal.  Thus  evaluating  the  correlation  without  chipping  codes  will 
not  give  us  the  clear  solution  we  are  seeking.  While  the  correlation  between  r{t)  and 
should  be  large  when  0,  matches  the  actual  received  direction  of  arrival  6inc,  there 
is  no  guarantee  that  for  any  other  arbitrary  angle  Q\,  the  correlation  would  not  be  large. 

More  generally,  we  can  clearly  see  the  advantage  of  the  chipped  sequence  by 
comparing  the  correlation  with  and  without  the  spreading  sequence.  In  the  non-spread 
approach,  the  correlation  sum  contains  cross  terms  as  seen  in  Equation  (42)  while  the 
SDMA  approach  produces  the  array  factor  in  which  the  cross  terms  are  zero  as  given  in 
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Equation  (36).  It  is  clear  that  the  orthogonal  spreading  codes  used  in  the  SDMA 
technique  are  integral  to  the  method. 

6.  Advantages  of  the  SDMA  Technique 

The  strength  of  this  technique  comes  from  the  ability  to  receive  several  signals 
and  “look”  in  almost  all  directions  nearly  simultaneously.  Because  the  chipping  of  the 
phase  of  the  received  signals  and  the  phase  offset  at  each  of  the  antenna  elements,  this 
essentially  makes  the  set  of  possible  antenna  patterns  independent  during  the  sampling 
period.  The  orthogonality  of  spreading  codes  between  antenna  elements  allows  the  actual 
signals  to  be  extracted  with  the  correct  direction  of  arrival.  The  signal  message  is 
extracted  from  the  phase  of  correlation  sequence.  Since  the  set  of  virtual  signals  is 
computed  once  and  stored  in  memory,  the  computational  cost  is  significantly  reduced  as 
compared  to  the  traditional  DO  A  methods,  such  as  conjugate  gradient  and  other  iterative 
matrix  solution  methods.  Additionally,  the  technique  provides  a  processing  gain  of 

Aff  [25]. 

Having  described  the  SDMA  method  we  now  demonstrate  through  simulation  the 
effectiveness  of  the  method  as  applied  to  a  uniform  linear  array  and  a  random  planar 
array  of  antenna  elements. 


D.  MATLAB  SIMULATION  RESULTS 

In  this  section  we  demonstrate  the  results  of  the  SDMA  technique  using 
MATLAB.  First  we  apply  the  technique  to  an  11  element  ULA  and  demonstrate  the 
importance  of  choosing  the  optimum  spreading  codes.  Having  selected  the  Walsh- 
Hadamard  codes  we  then  show  the  affect  of  the  number  of  chips  used  to  spread  the 
incoming  signal.  We  then  apply  the  technique  to  the  random  planar  array  and  simulate 
solving  for  a  DOA  and  then  extracting  a  phase  encoded  message  from  the  detected  signal. 
This  is  a  simulation  of  the  WSN  B  of  Figure  20  receiving  a  message  from  WSN  A  and 
then  decoding  the  location  message.  In  the  message  extraction  we  simulate  a  white 
Gaussian  noise  channel. 
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1. 


Uniform  Linear  Array 


In  this  simulation  we  use  the  SDMA  method  to  calculate  the  DOA  of  two 
simultaneous  test  signals.  The  ULA  contains  11  z-direction  oriented  dipole  elements 
spaced  one-half  wavelength  apart  and  is  used  to  detect  test  signals  arriving  from  two 
different  directions  -30°  and  +45°.  The  array  factor  for  this  ULA  is  given  simply  as  the 
sum  of  the  array  factors  from  each  incident  angle.  For  an  N  element  ULA,  the  array  factor 
for  a  single  incident  wave  from  direction  0inc  is  [18]: 

sin  ( Nkd  ( sin  0  -  sin  Qmc 
N  sin  ( kd  ( sin  0  -  sin  6inc  )  j 

and  the  complete  array  factor  for  two  incident  signals  is  the  sum  of  both  array  factors: 

AfT.Je)=AFx.(e)+AFt4S.(e) 

Again,  the  incident  “look”  angles  can  be  thought  of  as  the  beam  steering  angles.  Here, 
the  incident  beams  are  restricted  to  the  xz  plane.  The  geometry  used  in  this  simulation  is 
shown  in  Figure  25. 


Figure  25.  Geometry  for  ULA  with  Signals  Incident  from  -30°  and  +45° 

The  choice  of  spreading  codes  has  an  influence  on  the  quality  of  the  DOA 
determination.  In  Figure  26,  we  compare  the  result  of  the  correlation  using  two  different 
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types  of  spreading  codes,  the  Walsh-Hadamard  codes  and  a  pseudorandom  noise  code 
sequence.  The  correlation  magnitudes  and  the  array  factor  are  plotted  against  the  angle  of 
incidence  6. 


Figure  26.  SDMA  -  Comparison  of  PN  and  Walsh  Spreading  Codes  Using  32  Chips 

1 1  Element  ULA  with  Signals  Incident  from  -30°  and  +45°. 

The  green  line  is  the  array  factor.  In  this  comparison,  we  used  32  chips  for  both 
the  pn  sequence  and  the  Walsh-Hadamard  codes  and  4  samples  of  the  signal  per  chip.  A 
low  number  of  chips  was  chosen  in  order  to  make  the  difference  in  spreading  code  choice 
more  obvious.  The  blue  line  is  the  correlation  magnitude  calculated  using  the  Walsh- 
Hadamard  spreading  codes  and  the  magenta  line  is  the  correlation  magnitude  calculated 
using  the  pn  sequence.  In  this  figure,  Walsh-Hadamard  codes  clearly  work  better  than  a 
pn  sequence.  The  pn  codes  are  adequate  at  determining  the  DOA  of  the  signal,  but  the 
sidelobes  do  not  match  closely  with  those  of  the  theoretical  array  factor.  For  a  random 
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array,  where  the  sidelobes  do  not  decrease  as  strongly  as  the  ULA,  this  suggests  that  the 
Walsh-Hadamard  codes  are  a  better  choice. 

We  now  examine  the  effect  of  the  number  of  chips  used.  In  Figure  27,  we 
increased  the  number  of  chips  used  from  32  to  64  with  the  same  number  of  samples  per 
chip,  array  geometry  and  incident  signals  as  in  Figure  25.  The  array  factor  for  the  ULA 
is  plotted  in  green  and  the  magnitude  of  the  correlation  A,  is  plotted  in  blue.  We  see  from 
this  result  that  the  SDMA  technique  produces  a  correlation  in  excellent  agreement  with 
the  theoretically  produced  array  factor,  clearly  showing  the  DOA  of  the  two  incident 
signals  at  8  =  -30°  and  +45°  and  accurately  matching  the  mainlobe  as  well  as  all  the 
sidelobes. 


Figure  27.  DOA  Determination  of  Two  Incident  Signals  Using  SDMA  on  an  1 1  Element 
Array  with  One  Half  X  Spacing  and  64  Chip  Walsh-Hadamard  Code 


From  Figures  26  and  27,  we  see  that  the  choice  of  spreading  code  and  the  number 
of  chips  are  important.  In  comparison  to  Figure  26,  we  see  that  64  chips  produce  a  much 
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closer  match  to  the  theoretical  array  factor  than  32  chips  when  using  Walsh-Hadamard 
codes  on  the  same  array.  Additional  comparisons  (not  shown)  reveal  that  increasing  the 
number  of  chips  to  128  or  greater  had  only  marginal  increase  in  the  correlation  magnitude 
as  compared  to  theory.  It  appears  that  64  chips  are  adequate  for  DOA  determination. 

Having  demonstrated  that  the  SDMA  receiver  performance  agrees  with  theory  for 
the  uniformly  spaced  linear  array,  we  now  extend  the  technique  to  the  problem  of  the 
planar  array  of  randomly  placed  elements. 


2.  MATLAB  Simulation  Results  -  SDMA  and  Random  Planar  Array 


In  section  A  of  this  chapter,  we  discussed  the  need  for  the  receiving  array  to 
determine  both  the  DOA  of  the  searching  signal  and  the  range  to  the  searching  WSN. 
Here  we  apply  the  SDMA  technique  to  the  random  planar  array  of  dipoles  to  determine 
the  DOA  of  a  test  signal  and  compare  the  gain  pattern  from  the  SDMA  technique  to  that 
of  the  array  as  calculated  in  Chapter  III.  Having  found  the  DOA,  we  simulate  the 
decoding  of  a  sample  message  representing  bits  of  a  data  stream  from  the  searching 
WSN.  The  method  for  calculating  the  contribution  of  each  array  element  to  the 
correlation  value  for  a  random  planar  array  is  exactly  analogous  to  that  of  Chapter  III; 
instead  of  using  a  geometric  expression  for  the  location  of  a  receive  element  we  use  the 
exact  location  as  was  done  in  Equation  (14).  We  modify  Equation  (36)  for  the  random 
planar  array  as  follows: 


R.  =  oA)e 

n= 1  m= 1 


jtym  (6,  • A-  )eKO) 


t+T  .n  .n 

•  2  e  2  dt 


(46) 


where  the  element  phase  i//n(0o,0o)  is  now  determined  by  its  x-y  coordinate  and  the 


steering  angle  (0O.0O)  given  by: 


»//, (60,(j)0)=  (xn  sin 0 cos (j)  -  xn  sin 00  cos </>0 )  +  (y„  sin 9 cos 0  -  yn  sin 00  cos 0O ) 


a.  DOA  Determination 

We  now  apply  the  SDMA  technique  to  a  planar  array  of  randomly  placed 
sensors.  Figure  28a  depicts  an  overhead  view  of  the  elements  of  a  5/1  x  5  A  sensor  WSN 
containing  30  elements,  randomly  distributed  with  the  array  scanning  the  x-y  plane. 
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Figure  28b  depicts  the  simulated  test  signal  incident  from  9  =  90°  and  (f>  =  30°.  In  this 
diagram  R((9,</>)  is  shown  to  establish  the  coordinate  system  and  indicates  the  direction 
of  incidence  of  the  test  signal. 
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Figure  28.  Thirty  Element  Random  Planar  Array  Realization  and  Coordinate  System 


With  the  geometry  shown  in  Figure  28,  the  gain  pattern  of  the  random 
planar  array  was  calculated  using  Equations  (3)  and  (14)  and  compared  to  the  SDMA 
gain  pattern  calculated  using  Equations  (3)  and  (46).  Equation  (3)  is  the  element  pattern 
of  the  half  wave  dipole.  In  Figure  29,  the  green  line  shows  the  gain  pattern  of  the  random 
planar  array  with  the  incident  test  signal;  it  is  the  signal  the  SDMA  technique  must  locate. 
The  SDMA  gain  pattern  is  presented  in  blue  and  shows  excellent  agreement  with  the 
calculated  gain  pattern. 


65 


Figure  29.  Gain  Pattern  of  30  Element  Random  Planar  Array,  Steered  to  30°  with  SDMA 

Estimation  of  Incident  Signal  Pattern. 

As  shown  in  Chapter  III,  the  random  array  produces  a  main  lobe  in  the 
steered  direction  with  the  sidelobes  displaying  the  characteristic  random  nature  below  a 
pedestal.  The  30  elements  produce  a  peak  gain  of  17.8  dB,  the  beamwidth  is  8  degrees 
and  the  highest  side  lobe  rises  to  8.3  dB  below  the  peak.  These  values  are  in  agreement 
with  the  theoretical  results  presented  in  Chapter  III.  The  SDMA  pattern  accurately  shows 
the  main  lobe  at  30°  and  all  of  the  sidelobe  peaks.  The  angle  where  the  nulls  of  the  gain 
pattern  occur  are  accurately  shown  with  the  SDMA  technique,  but  the  absolute  values  do 
not  exactly  agree  with  the  theoretical  values. 

These  results  validate  the  SDMA  technique  as  a  method  for  blindly 
determining  the  DOA  when  applied  to  the  random  planar  array.  It  shows  excellent 
agreement  between  the  theoretical  array  patterns  and  those  produced  by  MATLAB.  This 
will  allow  the  receiving  WSN  to  almost  instantaneously  determine  the  DOA  of  the 
searching  WSN. 

Having  determined  the  DOA  of  the  searching  WSN,  the  task  of  the 
receiving  WSN  is  to  extract  the  location  information  from  the  message  signal  sent  by  the 
searching  WSN. 
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b.  Message  Extraction 

One  of  the  primary  goals  of  the  receiving  array  is  to  extract  from  the 
transmitted  signal  a  message  containing  the  location  of  the  searching  WSN.  In  order  to 
demonstrate  this,  we  will  simulate  the  baseband  message  signal  as  a  bit  sequence  and 
include  it  in  the  signal  incident  from  30°.  The  goal  is  to  transmit  an  8  bit  sequence  (Oil 
0  10  11)  mapped  to(l  -1  -1  1  -1  1  -1  -l)as  the  sample  message.  The  signal  phase  £(t) 
was  generated  with  additive  white  Gaussian  noise  (AWGN)  with  a  20  dB  signal  to  noise 
ratio  to  simulate  channel  interference.  Here  MATLAB  measures  the  power  in  £(t)  and 
adds  AWGN  at  20  dB.  The  signal  to  noise  ratio  is  the  ratio  of  the  power  in  the  signal  £  (t) 
to  the  power  in  the  AWGN. 

The  results  of  this  simulation  are  shown  in  Figure  30.  We  see  in  blue,  the 
instantaneous  baseband  signal  £(t)  with  AWGN.  In  green  is  the  phase  term  p,=30  of  the 

correlation  sequence  Rj  at  the  incidence  angle  0  =  30°.  The  time  sample  of  the  received 
signals  is  of  length  T  and  shown  in  red.  For  simplicity,  we  send  one  bit  per  time  period  T, 
although  this  need  not  be  the  case  [20,25].  We  see  that  the  signal  bit  sequence  can  be 
accurately  extracted  from  the  recovered  estimate  of  the  phase. 
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Figure  30.  SDMA  Message  Extraction  of  Sequence  [01101011] 
from  Noisy  Baseband  Signal,  SNR  =  20  dB 

In  order  to  test  the  strength  of  the  SDMA  technique,  we  increased  the  noise  level 
such  that  the  SNR  was  now  -5  dB.  Figure  30  shows  the  received  baseband  test  signal  in 
blue  encoding  the  same  bit  sequence  at  the  higher  noise  level.  The  recovered  signal  is 
shown  in  green.  Although  the  noise  level  appears  very  high,  the  SDMA  technique  is  still 
able  to  extract  the  data  in  the  incident  signal.  This  suggests  that  the  receiving  WSN  is 
capable  of  extracting  the  message  bits  under  high  noise  conditions. 
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Figure  31.  SDMA  Message  Extraction  of  Sequence  [01101011] 
from  Noisy  Baseband  Signal,  SNR  =  -5  Db 


The  direction  of  arrival  for  a  signal  from  a  searching  WSN  is  determined  by  the 
magnitude  of  the  correlation  sequence  and  the  message  itself  is  contained  within  the 
phase  of  the  correlation.  The  bit  stream  is  the  message  containing  the  “name”  and 
“location”  of  the  WSN  searching  for  the  adjacent  WSN.  Thus,  the  two  key  pieces  of 
information  required  to  establish  the  link,  the  bearing  and  range  can  now  be  determined. 
The  latitude  and  longitude  location  information  can  be  converted  into  a  range  using 
Equations  (18)  and  (19). 

The  results  shown  here  are  preliminary  to  demonstrate  the  effectiveness  of  the 
proposed  scheme.  Further  in  depth  simulation  and  comparison  to  gain  a  full 
understanding  and  scope  of  the  proposed  technique  is  required. 

In  this  chapter,  we  established  the  need  for  the  receiving  WSN  to  blindly 
determine  the  DOA  of  an  arriving  signal.  A  new  SDMA  receiver  was  described  and  its 
effectiveness  is  demonstrated  through  the  simulation  of  a  random  planar  array.  We 
showed  that  the  method  works  better  for  Walsh-Hadamard  spreading  codes  than  for  pn 
codes  generated  by  a  pseudorandom  number  generator  and  that  a  64  chip  sequence  is 
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sufficient  to  accurately  establish  the  DOA  and  extract  the  message.  Using  a  random 
planar  array,  we  validated  the  technique  by  determining  the  DOA  of  an  incident  test 
signal  and  then  extracting  a  message  data  stream  from  this  signal  in  a  noisy  channel  at 
SNRs  of  20  db  and  -5  dB.  From  the  DOA  and  the  location  information  within  the 
arriving  signal,  the  WSN  has  sufficient  knowledge  to  respond  to  the  query  of  the 
searching  WSN. 

Once  the  receiving  array  has  knowledge  it  is  being  queried  for  a  link,  the  task  is 
then  to  focus  its  array  towards  the  searching  array,  calculating  the  necessary  parameters 
required  for  the  link.  This  is  done  via  the  link  budget  analysis  as  was  described  in 
Chapter  II,  thus  completing  the  link  between  two  WSNs. 
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V.  CONCLUSIONS 


In  this  chapter,  we  present  a  summary  of  the  work  done,  highlight  the  significant 
results  and  suggest  topics  for  further  research. 


A.  SUMMARY  OF  RESEARCH 

In  this  thesis  we  explored  methods  for  forming  an  Over  The  Horizon  (OTH) 
communications  link  of  Wireless  Sensor  Networks  (WSNs)  by  enabling  each  WSN  to  act 
as  a  smart  antenna  array.  Methods  for  establishing  the  OTH  communications  link  via 
beam  forming  and  direct  sequence  spread  spectrum  Space  Division  Multiple  Access 
(SDMA)  were  presented  and  modeled  using  MATLAB.  Methods  for  forming  a  search 
beam  for  a  transmitting  WSN  and  determining  the  direction  of  arrival  of  the  received 
beam  for  a  receiving  WSN  were  presented. 

We  examined  the  searching  WSN  in  which  the  individual  sensor  nodes  operate  in 
concert  as  a  smart  antenna  array.  The  WSN  then  forms  and  transmits  a  search  beam 
around  the  horizon  in  order  to  illuminate  an  adjacent  WSN  and  establish  a  link.  Each 
WSN  was  modeled  as  a  uniformly  distributed  random  array  of  sensors  nodes.  Gain  and 
beamwidth  of  the  random  planar  arrays  were  calculated  and  validated  using  MATLAB 
simulation  and  compared  to  theoretical  results.  Through  MATLAB  simulation,  the  half 
power  beamwidth  (HPBW)  of  a  random  planar  array  as  a  function  of  WSN  size  and  the 
gain  as  a  function  of  the  number  of  elements  for  a  fixed  size  array  were  studied.  A  40 
element  random  planar  array,  comparable  to  a  36  element  periodic  planar  array  of  similar 
size,  was  simulated  to  demonstrate  the  search  capability.  . 

We  examined  the  proposed  communications  link  from  the  perspective  of  the 

receiving  array.  This  array  needs  to  determine  the  direction  of  arrival  (DOA)  and 

location  of  the  searching  WSN  in  order  to  form  an  optimum  link  without  a  reference  or 

pilot  signal.  A  SDMA  receiver  capable  of  detennining  the  DOA  and  location  of  the 

searching  WSN  was  described  and  its  effectiveness  demonstrated  through  the  simulation 

of  a  random  planar  array.  Using  a  random  planar  array,  we  validated  the  SDMA  method 
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by  determining  the  DOA  of  an  incident  test  signal  and  then  extracting  a  message  data 
stream  from  this  signal  in  a  noisy  channel.  From  the  DOA  and  the  location  information 
within  the  arriving  signal,  the  WSN  has  sufficient  knowledge  to  respond  to  the  query  of 
the  searching  WSN. 


B.  SIGNIFICANT  RESULTS 

We  investigated  the  random  planar  array  performance  in  terms  of  the  half  power 
beam  width,  the  gain  of  an  array  of  N  elements  and  side  lobe  performance.  The  objective 
was  to  apply  these  results  to  wireless  sensor  networks  to  first  form  a  search  beam  to 
locate  adjacent  WSNs,  then  form  a  narrow  high  gain  beam  to  form  a  communications 
link.  The  particular  WSNs  are  assumed  to  be  fixed  in  location,  but  the  inter-element 
spacing  follows  a  uniform  random  distribution  within  a  defined  boundary.  We 
demonstrated  that  we  could  select  a  wide  beam  for  the  search  pattern  and  simultaneously 
select  the  gain.  To  achieve  a  given  beamwidth,  we  can  choose  a  physical  size,  likely  a 
sub-section  of  the  WSN  under  consideration.  Within  this  physical  area,  we  have  control 
over  the  gain  by  choosing  the  number  of  elements  employed  in  the  array.  Alternately,  for 
a  given  required  gain  level  we  can  thin  the  array  (or  reduce  density)  without  significant 
loss  of  gain  or  the  effects  of  grating  lobes.  It  does  not  matter  which  particular  sensor 
nodes  within  the  WSN  are  used  as  long  as,  on  average,  they  are  randomly  distributed 
within  the  physical  area  chosen.  This  is  important  if  the  objective  is  to  distribute  the 
burden  of  transmitting,  and  thus  energy  consumption  among  the  nodes  within  the  WSN. 

We  established  the  need  for  the  receiving  WSN  to  blindly  determine  the  direction 
of  arrival  (DOA)  of  an  arriving  signal.  From  the  DOA  and  the  range  information 
encoded  within  the  arriving  signal,  the  WSN  has  sufficient  knowledge  to  build  its  reply  to 
the  query  of  the  searching  WSN.  Using  a  random  planar  array,  we  utilized  a  spread 
spectrum  base  SDMA  technique  to  determine  the  DOA  of  an  incident  test  signal  and 
simulated  extracting  a  message  data  stream  from  this  signal  in  a  noisy  channel  with  SNRs 
of  20  db  and  -5  dB,  accurately  extracting  the  data. 
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C.  TOPICS  FOR  FURTHER  RESEARCH 

In  the  process  of  carrying  out  this  work  many  issues  for  further  research  came  to 
light.  This  section  highlights  some  of  these  topics. 

In  this  work,  we  demonstrated  a  SDMA  method  to  determine  the  DOA  of  a  signal 
in  two  dimensions  assuming  vertical  polarization.  Indeed,  the  SDMA  method  shows 
promise  for  much  greater  capabilities.  In  this  thesis,  the  virtual  signals  were  only  used  to 
determine  the  DOA  across  the  horizon  but  could  be  extended  to  include  the  hemisphere. 
In  this  manner,  the  graphically  depicted  UAV  of  Figure  4  could  quickly  be  located,  or  in 
the  reciprocal  instance,  the  UAV  could  quickly  find  WSNs.  Once  located,  the  look 
directions  or  virtual  signals  used  for  tracking  quickly  collapses  to  just  a  few  vectors.  This 
type  of  three  dimensional  pattern  covering  the  hemisphere  above  the  array  is  often 
referred  to  in  the  literature  as  a  “pin  cushion”  pattern.  The  SDMA  technique  can  also 
search  for  multiple  or  elliptical  polarizations  by  creating  and  summing  the  orthogonal 
polarization  virtual  signals  [25]. 

In  this  work,  the  spreading  was  applied  only  to  the  receiving  WSN  to  determine 
the  DOA  and  extract  the  message.  Alternately,  the  searching  WSN  could  spread  its 
transmitted  search  signal  in  the  manner  described  in  Chapter  IV  on  SDMA  for  the 
receiving  array  leading  to  further  processing  gain. 

For  the  SDMA  receiver,  we  assumed  that  only  the  direct  signal  reached  the 
receiver.  For  a  terrestrial  communications  link,  there  would  certainly  be  a  multipath 
signal  fading  [16].  Future  work  may  investigate  the  receiver’s  capability  to  reject 
multipath  signals  and  increase  the  communication  link  range  or  bit  rate. 

In  this  work,  we  assumed  that  the  location  of  each  sensor  and  that  the  location 
error  could  be  estimated.  Appendix  A  addresses  the  methods  of  self-localization  of  the 
array,  however  the  impact  of  these  errors  to  the  SDMA  technique  were  not  addressed  in 
this  thesis.  Future  work  may  determine  the  impact  of  localization  errors  on  the  SDMA 
technique.  Additionally,  future  work  may  address  methods  for  combining  the  signals 
received  at  each  of  these  individual  sensor  node  locations. 
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The  path  loss  model  in  this  work  used  a  wireless  LAN  model  based  on  antenna 
heights  of  1  m  or  above.  A  better  path  loss  model  would  give  a  more  accurate  estimate  of 
the  margin  required  to  form  the  communications  link  between  two  adjacent  WSNs. 
Future  work  may  seek  to  develop  path  loss  models  for  ground  level  transmitter-receiver 
pairs  based  on  empirical  data  using  actual  sensor  motes  under  a  variety  of  operating 
conditions.  Applications  for  linking  together  WSNs  are  only  limited  by  the  imagination. 
The  WSNs  need  not  be  random  or  fixed  in  place  on  the  ground.  For  example,  one  WSN 
may  be  a  periodic  array  placed  on  the  side  of  a  building  and  the  other  could  be  an  array 
placed  on  the  top  or  side  of  a  truck.  As  the  truck  passes  a  station,  the  fixed  station  would 
search  for  the  truck’s  array,  establish  a  communications  link  and  pass  information 
wirelessly  to  the  fixed  node.  Since  the  link  is  established  autonomously,  the  truck 
doesn’t  need  to  stop,  or  indeed  perform  any  action.  A  future  effort  may  examine  such 
applications  that  have  potential  military  value. 
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APPENDIX  A.  WIRELESS  SENSOR  NODE  LOCALIZATION 


A  basic  requirement  of  creating  a  smart  antenna  from  a  WSN  is  the  immediate 
task  of  forming  a  local  network  of  these  nodes,  identifying  each  other’s  location  and 
relative  geometry.  This  map  of  the  geometry  is  created  during  the  self  organization  of  the 
WSN.  Once  the  local  geometry  of  the  WSN  is  determined,  it  will  then  be  possible  to 
form  a  beam  collectively  and  search  for  other  WSNs.  This  section  deals  with 
determining  the  local  wireless  sensor  network  geometry. 

It  is  assumed  that  the  sensor  nodes  are  deployed  or  dispersed  in  a  random  fashion 
in  the  area  of  interest  and  that  all  are  properly  oriented,  i.e.  they  are  not  laying  on  their 
side  or  upside  down,  but  have  the  transmitting  antenna  pointing  vertically.  For  the 
purposes  of  this  work,  the  sensors  in  a  network  are  assumed  to  be  randomly  distributed 
with  a  uniform  probability  density  function  (pdf). 


A.  LOCALIZATION  METHODS 

So  important  is  the  need  to  determine  local  geometries  of  an  ad  hoc  wireless 
network,  extensive  research  has  been  conducted  on  this  topic  and  the  corresponding 
literature  covers  a  multiplicity  of  techniques  [4,30,31,32,33].  Here  we  present  a  brief 
description  of  the  methods  for  determining  the  geometry  of  an  ad  hoc  wireless  sensor 
network.  Due  to  the  available  methodologies  and  advancing  technology,  we  assume  in 
this  thesis  that  each  local  WSN  is  self  aware  of  its  absolute  location  and  orientation  as 
well  as  the  geometry  of  the  individual  nodes  that  compose  the  network. 

Self  organization  allows  the  WSN  to  determine  its  own  geometry.  These 
techniques  typically  estimate  distance  and  angle  between  nodes,  communicate  them  to 
the  neighboring  nodes,  then  adjust  or  rotate  individual  coordinate  systems  to  obtain  a 
single  global  coordinate  system.  Magnani  and  Leung  [4]  discuss  one  such  algorithm. 

Localization  methods  generally  fall  into  4  major  categories  and  are  described 

here. 
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1. 


GPS  Localization 


All  the  sensor  nodes  have  on-board  GPS  receivers  and  simply  communicate  this 
to  adjacent  nodes,  forming  the  table  of  locations  and  or  separation  distances. 

2.  Signal  Techniques 

These  techniques  rely  on  nodes  estimating  the  location  of  other  nodes  by 
measuring  a  signal  parameter  [30],  One  such  method  measures  the  received  signal 
strength  as  compared  to  a  known  transmitted  value  called  Received  Signal  Strength 
Indicator  (RSSI)  or  Relative  Signal  Strength  (RSS).  Other  techniques  are  Time  of  Arrival 
(ToA),  Time  Difference  of  Arrival  (TDoA)  and  Angle  of  Arrival  (AoA)  estimation.  For 
the  AoA  case,  the  sensor  node  must  be  equipped  with  an  antenna  that  is  able  to  measure 
direction. 


3.  Anchor  Nodes 

Anchor  nodes  or  “base  stations”  of  known  position  are  distributed  across  a  sensor 
network  [30,3 1].  Other  nodes’  positions  are  triangulated  based  on  the  known  locations  of 
the  anchors.  Careful  selection  of  the  base  stations’  position  distribution  and  precise 
knowledge  of  their  location  are  crucial  to  minimizing  position  error  propagation  as 
subsequent  node  locations  are  calculated  [31].  Additionally,  the  triangulation  error  in 
locating  a  particular  node  can  be  reduced  by  dispersing  the  base  stations  at  the  perimeter 
edges  of  the  sensor  field  as  opposed  to  picking  3  base  stations  that  are  close  together. 

4.  GPS  Free  Methods 

Numerous  references  are  available  for  the  localization  of  WSN  in  the  absence  of 
GPS  data  [4,32,33].  These  methods  require  a  signal  technique  and  the  determination  of 
local  and  absolute  coordinate  systems.  Reference  [31]  shows  that  multi-path  errors 
dominate  the  measurement  noise  and  subsequently  tend  to  be  the  main  cause  of  error  in 
range  estimation. 

Ji  and  Zha  [30]  discuss  the  method  of  multidimensional  scaling  as  technique  to 
mitigate  the  problems  of  cumulative  measurement  error.  This  technique  builds  a  series  of 
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local  maps  of  adjacent  sensors  along  a  route  from  a  starting  anchor  to  an  ending  anchor. 
This  method  may  be  of  use  for  the  larger  problem,  discussed  in  Chapter  II  on  building  the 
OTH  network  of  WSNs. 

In  this  work,  we  assume  the  following  about  the  WSNs:  the  WSN  is  location 
aware  with  respect  to  its  absolute  location  and  orientation;  the  relative  positions  of  each 
of  the  sensor  nodes  comprising  the  network  are  known;  the  sensor  nodes  are  uniformly 
distributed  within  a  definable  boundary;  and  the  sensor  nodes  are  able  to  synchronize 
internal  clocks  either  to  a  master  node  or  to  a  GPS  reference. 


77 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


78 


APPENDIX  B.  MATLAB  CODE 


This  appendix  contains  selected  MATLAB  code  used  to  generate  the  figures  in  this 
thesis. 

A.  Comparison  of  Power  Pattern  of  Linear  Array  and  Average  Power  Pattern 
of  the  Same  Size  Random  Array 

This  code  was  used  to  generate  Figure  12. 


%  Chris  Taylor  Modification  to  Array  factor  code  presented  as  a  solution  to 
%  Homework  4.1  of  below  textbook. 

%  This  code  produces  both  the  Array  factor  for  a  uniform  linear  array  and 
%  the  envelope  from  a  random  array  of  the  same  size. 

%  This  code  also  used  in  the  Figure  9  on  grating  lobes. 

%  Feb  08  Naval  Postgraduate  School,  Monterey  CA 

% 

%  Algorithm  for  ULA  by  Frank  B.  Gross  "Smart  Antennas  for  Wireless 
%  Communications",  c  2005,  Chapter  4,  Homework  solutions. 

% 

%  Matlab  m-file  for  Problem  4.1 

d  =  [.5  0.5  0.5];  %  d/lambda 

%N  =  [888];  theta  =  -pi/2: .01: pi/2; 

N  =  [5  8  20];  theta  =  -pi:. 01: pi; 

% - Create  Normalized  Array  Factors - % 

AFnl  =  l/N(l)*sin(N(l)*pi*d(l)*sin(theta)) . /sin(pi*d(l)*sin(theta)) ; 

AFn2  =  l/N(2)*sin(N(2)*pi*d(2)*sin(theta)) . /sin(pi*d(2)*sin(theta)) ; 

AFn3  =  l/N(3)*sin(N(3)*pi*d(3)*sin(theta)) . /sin(pi*d(3)*sin(theta)) ; 

% - linear  Array  Power  Pattern - 

PAFnl  =  AFnl. *conj (AFnl);  %PAFnl=10*logl0(abs(PAFnl)) ; 

PAFn2  =  AFn2.*conj(AFn2);  %PAFn2=10*logl0(abs(PAFn2)) ; 

PAFn3  =  AFn3.*conj(AFn3);  %PAFn3=10*logl0(abs(PAFn3)) ; 

% — Random  Array  Predicted  Power  Pattern  Envelope - % 

RAFN1  =  (AFnl.*conj(AFnl)).*(l-(l/N(l)))+  1/N(1)  ; 

RAFN2  =  (AFn2.*conj(AFn2)).*(l-(l/N(2)))+  1/N(2)  ; 

RAFN3  =  (AFn3.*conj(AFn3)).*(l-(l/N(3)))+  1/N(3)  ; 

thetad=theta*180/pi ; 

% - Plot  Results - %  4-1 

subplot(311) ,  plot(thetad, 10*logl0(PAFnl) , ' b' , thetad, 10*logl0(RAFNl) , ' -c ' ) ,  grid  on 
title([  {'N  Element  linear  Array  Power  Pattern'}> • • • 

{'  vs  Average  Power  Pattern  of  N  Element  Random  Array'}],... 

' FontSize' ,14, ' FontWeight' , 'bold') 
ylabel(  \fontsize{14}AFA2(\theta)  dB ' ) 
xlabel(  \fontsize{14}\theta  (deg)  ) 

legend('N  =  5,  d/\lambda  =  0.5',  'Random  Array  Ave  PP'), 
axis([0  90  -20  0]) 

subplot(312) ,  plot(thetad,10*logl0(PAFn2), ' r' , thetad, 10*logl0(RAFN2), ' -m'),  grid  on 
ylabel(  \fontsize{14}AFA2(\theta)  dB ' ) 
xlabel(  \fontsize{14}\theta  (deg)  ) 
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legendC'N  =  8,  d/\lambda  =  0.5',  'Random  Array  Ave  PP'), 
axis([0  90  -30  0]) 

subplot(313) ,  plot(thetad,10*logl0(PAFn3), ' b ' , . . . 

thetad,10*logl0(RAFN3), ' g ' ),  grid  on 

ylabel('\fontsize{14}AFA2(\theta)  dB ' ) 


xlabel( ' \fontsize{14}\theta  (deg) ' ) 

legendC'N  =  20,  d/\lambda  =  0.5',  'Random  Array  Ave  PP'), 
axis([0  90  -40  0]) 


B.  Random  Planar  Array  Half  Power  Beam  Width  vs  Array  Size  and 

Random  Planar  Array:  Maximum  Gain  vs  Array  Size,  Const  Number  of 
Elements 

This  code  was  used  to  generate  Figures  14  and  16 

%  This  code  computes  far-field  gain  pattern 
%  for  planar  array  composed  of  N  vertical  half-wave 
%  dipoles  at  wavelength  normalized  coordinates  [x(n),y(n)] 

%  with  complex  I(n)  currents. 

% 

%  Original  code  for  period  arrays  by  M.  A.  Morgan  29  Oct  07 

% 

%  Modified  by  C.E.  Taylor  Dec  07  to  use  Random  Planar  Array 
%  **  This  version  calculates  half  power  BW  vs  array  size  ** 

% 

%  This  version  uses  Uniformly  current  amplitudes  (In=l  for  all  elem) 

% 

clear  all 
elf  reset 


%  Setting  up  3-D  pattern  grid 

Ntheta=91;  Nphi=361; 

theta=linspace(0 , pi/2 , Ntheta) ;  phi=linspace(-pi , pi , Nphi) ; 

[Phi , Theta] =meshgrid(phi , theta); 

CT=cos(Theta) ;  ST=sin(Theta) ;  CP=cos(Phi);  SP=sin(Phi); 

Rx=ST.*CP;  Ry=ST.*SP;  Rz=CT ;  %  unit  radial  vector  components 

%  Flalf-wave  dipole  pattern  for  Iin=lA  (eps  handles  theta=0  indeterminent) 
E0  =  60*cos(pi*CT./2)./(ST+eps); 


HPBW_ne  =  zeros(l,50); 
Gain_ne  =  zeros(l,50); 
PeakSLL_ne  =  zeros(l,50); 
AveSLL_ne  =  zeros(l,50); 


%  Flalf  Power  Beam  Width  vs  array  size 
%  Gain  vs  array  size 
%  Peak  SLL  vs  size 
%  Ave  SLL  vs  array  size 


% - Build  Random  Array - 

lambda  =  1.0; 
k  =  2*pi/lambda; 

%array_size  =  5.0*lambda;  %  in  wavelengths  (lambda)  square  grid 


nruns  =  35; 
for  run=l: nruns; 


fprintf(  :  run:  %2.0f\n',  run); 
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sp  so  so  so  sp 

OS  OS  OS  OS  OS 


for  array_size  =  (1: 20)*lambda; 
num_elements  =  20; 

Nx  =  num_elements ; 

Ny  =  num_elements ; 
sepn_min  =  0; 

xp  =  array_size  *  rand(l,num_elements); 
yp  =  array_size  *  rand(l,num_elements); 

% - Stats  of  Random  Array - 

min_spacing  =  0.1*lambda;  %  l/10th  wavelength 
while  (sepn_min<  min_spacing)  %  min  spacing  =  .1  wavelength 

rp  =  [xp;yp] ' ; 

sepn  =  pdist(rp);  %  dist  between  all  elements  in  Random  Array 
sepn_min  =  min(sepn);  %  find  minimum  separation 

zsepn  =  squareform(sepn) ;  %  convert  linear  sepn  array  to  square 

[Rmin  Cmin]  =  find(zsepn==sepn_min,l,  'first');  %returns  index 
fprintf(['Min  separation  is  %6.3f  wavelengths,'... 

'between  elements  %2.0f  and  %2.0f  \n '],... 
sepn_min,  Rmin,  Cmin) 

if  (sepn_min<  min_spacing)  %  then  Move  one  of  them 
xp(Rmin)=array_size  *  rand; 
yp(Rmin)=array_size  *  rand; 

end 
end 

fprintf(  'Element  %2.0f  located  at  %8.3f,  %8.3f\n' ,Rmin,rp(Rmin, :)  ) 
fprintf(  'Element  %2.0f  located  at  %8.3f  %8.3f\n' ,Cmin, rp(Cmin, :)  ) 

10=1;  %  assuming  uniform  current  amplitudes 

%phid0=input( ' Enter  peak  gain  azimuth  angle  in  degrees  from  x-axis:  ') 

phid0  =  0; 

phi0=pi*phid0/180;  %convert  to  radians 

%  Phasing  I(Nx,Ny)  array  elements  to  point  array  in  phi0  direction 
Ip=  10  *  exp(- j*k*(cos(phi0)*xp+sin(phi0)*yp)) ; 

%  Summing  far-field  pattern 

F=zeros(size(E0)) ; 
for  n=l:num_elements 
xpRx  =  xp(n)*Rx; 
ypRy  =  yp(n)*Ry; 

F  =  F  +  Ip(n)*  exp( j*k*(xpRx+ypRy)) ; 

end 

%  Normalized  power  pattern 
F2  =  abs((E0.*F)/(num_elements)).A2; 

DTheta=theta(2)-theta(l) ;  DPhi=phi(2)-phi(l) ; 

IntF2=DTheta*DPhi*sum(sum((F2) . *sin(Theta))) ; 

Gain=4*pi*F2/IntF2; 

Gmax=max(max(Gain)) ;  GmaxdB=10*logl0(Gmax) ; 

%  Horizontal  pattern  (the  x-y  plane)  and  3dB  beamwidth 
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n=findfphi==0) ; 
m=findftheta==pi/2) ; 
G3dB=Gmax/2; 

GH=Gain(m, : ) ; 

mf=findfGH  >=  G3dB,l, 'first  ); 
mp=findfGH  >=  G3dB,l, 'last'); 
BWH=180*Cphifmp)-phifmf))/pi ; 
GHdB=10*logl0(GH); 


% -  calc  Side  Lobe  levels  - 

BW  =  imregionalmaxfGHdB) ;  %MATLAB  rtn  that  finds  local  maxima->BW  array  "l"s 
A1  =  findfBW);  %index  pointers  to  local  Maxima 

[bb,ix]=sortfGHdBffindfBW)), 'descend');  %  sort  by  highest  peak->lowest 
SLL=  bbfl)-bb;  %  Side  Lobe  levels,  down  from  peak 

SLi=  Alfix);  %  index  to  peaks 

SLa=  180*phi(SLi)/pi ;  %  Angle  where  peaks  occur 

%  fprintff ' First  Sidelobe  Level:  %6.2f  below  max  at  phi=  %4.0f  \n' ,SLLf2),SLaf2)) 

%  fprintfC ' Second  Sidelobe  Level:  %6.2f  below  max  at  phi=  %4.0f  \n' ,SLL(3),SLa(3)) 

%  fprintf('Mean  Sidelobe  Level:  %6.2f  below  max  \n' ,meanCSLLC2:length(Al)))) 

HPBW_ne(array_size)  =  HPBW_ne(array_size)+BWH  ;  %  Gain  vs  number  elements  array 

theta3db  (array_size)  =  50.76  /  array_size;  %  Theory 

Gain_ne(array_size)=  Gain_ne(array_size)  +  GmaxdB; 

PeakSll_ne(array_size)=  PeakSLL_ne(array_size)+SLL(2) ;  %  Peak  SLL  vs  num  elements 

AveSLL_ne(array_size)=AveSLL_ne(array_size)+mean(StL(2:length(Al)));  %  Ave  SLL 

vs  num  elements 

end  %  array_size 


end  %  runs 

plot(l:array_size,  HPBW_ne(l:array_size)/nruns,  '  -A' ,  l:array_size,theta3db) 

tl=['Random  Planar  Array  num2str(num_elements) , '  Elements']; 

t2='HPBW  vs  Array  Size'; 

title({tl,t2}); 

xlabelf  Array  Size  \lambda') 

ylabelf  Half  Power  Beam  Width  Ao') 

grid  off 

legendf  Random  Array' , 'Analytic' ) 

%axis([  1  30  0  360]); 

figure; 

plot(l:array_size,  Gain_ne(l:array_size)/nruns, . . . 

l:array_size,  Gain_ne(l:array_size)/nruns  -  PeakSlL_ne(l:array_size)/nruns, . . . 
l:array_size,  Gain_ne(l:array_size)/nruns  -  AveSLL_ne(l:array_size)/nruns, . . . 
l:array_size,  Gain_ne(l:array_size)/nruns-13)  %  13  =  -101ogl0(20  elements) 
tl=['Random  Planar  Array  ',  num2str(num_elements) , '  Elements']; 
t2='Maximum  Gain  vs  Array  Size'; 
title({tl,t2}); 
axisC  [1  20  0  23]); 
grid  off 

xlabelf  Array  Size  \lambda'); 
ylabelf  Power  f dB) ' ) ; 

legendf 'Max  Gain',  'Peak  SLL',  'Ave  SLL', 'Ave  SLL-Theory' , 'Location' , 'Northwest') 
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C.  Random  Planar  Array:  Maximum  Gain  vs  Number  of  Elements 

This  code  was  used  in  Figure  15. 

%  This  code  computes  far-field  gain  pattern 
%  for  planar  array  composed  of  N  vertical  half-wave 
%  dipoles  at  wavelength  normalized  coordinates  [x(n),y(n)] 

%  with  complex  I(n)  currents. 

% 

%  By  M.  A.  Morgan  29  Oct  07 

% 

%  Modified  by  C.E.  Taylor  Dec  07  to  use  Random  Planar  Array 
%  **  This  version  calculates  gain  vs  number  of  elements  ** 

% 

%  This  version  uses  Uniformly  current  amplitudes  (In=l  for  all  elem) 

% 

clear  all 
elf  reset 

%  Setting  up  3-D  pattern  grid 

Ntheta=91;  Nphi=361; 

theta=linspace(0 , pi/2 , Ntheta) ;  phi=linspace(-pi , pi , Nphi) ; 

[Phi , Theta] =meshgrid(phi , theta); 

CT=cos(Theta) ;  ST=sin(Theta) ;  CP=cos(Phi);  SP=sin(Phi); 

Rx=ST.*CP;  Ry=ST.*SP;  Rz=CT ;  %  unit  radial  vector  components 

%  Half-wave  dipole  pattern  for  Iin=lA  (eps  handles  theta=0  indeterminent) 

E0  =  60*cos(pi*CT./2)./(ST+eps); 

Gain_ne  =  zeros(l,50);  %  Gain  vs  number  elements  array 

PeakSLL_ne  =  zeros(l,50);  %  Peak  SLL  vs  num  elements 
AveSLL_ne  =  zeros(l,50);  %  Ave  SLL  vs  num  elements 

% - Build  Random  Array - 

lambda  =  1.0; 
k  =  2*pi/lambda; 

array_size  =  5.0*lambda;  %  in  wavelengths  (lambda)  square  grid 


nruns  =  75; 
for  run=l: nruns; 

fprintf(  '  run:  %2.0f\n',  run); 

for  num_elements  =  5:1:50; 

Nx  =  num_elements ; 

Ny  =  num_elements ; 
sepn_min  =  0; 

xp  =  array_size  *  rand(l,num_elements); 
yp  =  array_size  *  rand(l,num_elements); 

% - Stats  of  Random  Array - 

min_spacing  =  0.1*lambda;  %  l/10th  wavelength 


while  (sepn_min<  min_spacing)  %  min  spacing  =  .1  wavelength 
rp  =  [xp;yp] ' ; 

sepn  =  pdist(rp);  %  dist  between  all  elements  in  Random  Array 
sepn_min  =  min(sepn);  %  find  minimum  separation 

zsepn  =  squareform(sepn) ;  %  convert  linear  sepn  array  to  square 

[Rmin  Cmin]  =  find(zsepn==sepn_min,l,  'first');  %returns  index 
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\o  \o  \o  \p  \o  so  so  so 

OS  OS  OS  OS  os  os  o\  OS 


fprintf(['Min  separation  is  %6.3f  wavelengths,'... 

'between  elements  %2.0f  and  %2.0f  \n '],... 
sepn_min,  Rmin,  Cmin) 

if  (sepn_min<  min_spacing)  %  then  Move  one  of  them 
xp(Rmin)=array_size  *  rand; 
yp(Rmin)=array_size  *  rand; 

end 

end 

fprintf(  'Element  %2.0f  located  at  %8.3f,  %8.3f\n' , Rmin, rp(Rmin, :)  ) 
fprintf(  'Element  %2.0f  located  at  %8.3f  %8.3f\n' ,Cmin,rp(Cmin, :)  ) 

10=1;  %  assuming  uniform  current  amplitudes 

%phid0=input( ' Enter  peak  gain  azimuth  angle  in  degrees  from  x-axis:  '); 

phid0  =  0; 

phi0=pi*phid0/180;  %convert  to  radians 

%  Phasing  I(Nx,Ny)  array  elements  to  point  array  in  phi0  direction 
Ip=  10  *  exp(- j*k*(cos(phi0)*xp+sin(phi0)*yp)) ; 

%  Summing  far-field  pattern 

F=zeros(size(E0)) ; 
for  n=l : num_elements 
xpRx  =  xp(n)*Rx; 
ypRy  =  yp(n)*Ry ; 

F  =  F  +  Ip(n)*  exp( j*k*(xpRx+ypRy)) ; 

end 

%  Normalized  power  pattern 

F2  =  abs((E0.*F)/(num_elements)).A2; 

DTheta=theta(2)-theta(l) ;  DPhi=phi(2)-phi(l) ; 

IntF2=DTheta*DPhi*sum(sum(CF2) . *sin(Theta))) ; 

Gain=4*pi*F2/IntF2; 

Gmax=max(max(Gain)) ;  GmaxdB=10*logl0(Gmax) ; 

%  Horizontal  pattern  (the  x-y  plane)  and  3dB  beamwidth 

n=find(phi==0) ; 
m=find(theta==pi/2) ; 

G3dB=Gmax/2; 

GH=Gain(m, : ) ; 

mf=find(GH  >=  G3dB,l, 'first'); 
mp=find(GH  >=  G3dB,l, 'last'); 

BWH=180*(phi(mp)-phi(mf))/pi ; 

GHdB=10*logl0(GH) ; 


% -  calc  Side  Lobe  levels  - 

BW  =  imregionalmax(GHdB) ;  %MATLAB  rtn  that  finds  local  maxima->BW  array  "l"s 

A1  =  find(BW);  %index  pointers  to  local  Maxima 

[bb,ix]=sort(GHdB(find(BW)), 'descend');  %  sort  by  highest  peak->lowest 

SLL=  bb(l)-bb;  %  Side  Lobe  levels,  down  from  peak 

SLi=  Al(ix);  %  index  to  peaks 

SLa=  180*phi(SLi)/pi ;  %  Angle  where  peaks  occur 

fprintf( ' First  Sidelobe  Level:  %6.2f  below  max  at  phi=  %4.0f  \n' ,SLL(2),SLaC2)) 
fprintfC Second  Sidelobe  Level:  %6.2f  below  max  at  phi=  %4.0f  \n' ,SLL(3),SLa(3)) 
fprintf( 'Mean  Sidelobe  Level:  %6.2f  below  max  \n' ,mean(SLL(2 : length(Al)))) 
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Gain_ne(num_elements)=  Gain_ne(num_elements)  +  GmaxdB; 

PeakSLL_ne(num_elements)=  PeakSLL_ne(num_elements)+SLL(2) ;  %  Peak  SLL  vs  num 

elements 

AveSLL_ne(num_elements)=AveSLL_ne(num_elements)+mean(SLL(2: Length(AL))) ;  %  Ave 

SLL  vs  num  elements 


end  %  num_elements 


end  %  runs 


ne  =  l:num_elements; 
Gain_ne=Gain_ne/nruns ; 

PSSL  =  Gain_ne-PeakSLL_ne/nruns ; 

ASSL  =  Gain_ne-AveSLL_ne/nruns ;  % 

Gth=  Gain_ne+  10*Logl0(l./ne);  % 

Gain_th2  =  10*Logl0(ne)+10*Logl0(l . 6) ; 


Peak  SLL  actual 

Ave  SLL  actual 

Theory  ASLL  (Steinberg) 

%  Theory:  Gain  of  N  dipoles,  incl  gain  of  dipole 


plot(ne,  Gain_ne,  ne,PSSL,... 

ne,ASSL,  ne,Gth,  ne,  Gain_th2  ) 

tl=['Random  Planar  Array  ',  num2str(array_size) , '  \lambda  square']; 
t2='Maximum  Gain  vs  Number  of  Elements'; 
title({tl,t2}); 
axis(  [5  50  0  23]); 

xlabel(  Number  of  Elements  in  Array'); 
ylabel(  Power  (dB)'); 

legend('Gain' ,  'Peak  SLL',  'Ave  SLL',  'Theory-Ave  SLL',... 

'Theory-N  dipoles',  'Location' , 'Northwest') 
grid  off 


D.  DOA  determination  for  the  Uniform  Linear  Array,  with  two  incident  signals. 

This  code  used  in  Figure  23,  as  an  example  of  the  SDMA  technique. 


%  Modified  by  Chris  Taylor  Jan  2008  to  include  pn  seq  generation  or  Walsh 
%  codes,  samples  per  chip,  and  two  incident  beams 

%  This  is  a  simulation  to  model  the  Elam  SDMA  technique 
%  Developed  by  Frank  B.  Gross,  November  11,  2004 

%  define  the  sample  baseband  modulation  initially  as  a  1kHz  phase  modulation  tone 
%  WALSH  FUNCTIONS  WORK  BETTER  THAN  PN  SEQUENCES , Can  also  use  Gold  codes 
clear;  clc; 

f=lE3; 

T=l/(4*f); 

Nchips=64; 

SPC  =4;  %  Samples  Per  Chip  (4  min) 

tend=SPC*Nchips ;  %  time  snapshot  end 

t=[0 : (tend-1)] *T/(tend-l) ; 
zeta=sin(2*pi*f*t) ;  %  modulation  m 

%m=zeros(l , Length(m)) ; 

%  plot  array  factor  for  linear  array 

N=ll;  %  number  of  receive  elements  to  give  10  degree  beamwidth 
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d=.5;  %  spacing  in  wavelengths  between  adjacent  array  elements 


%  generate  pseudo-random  chip  codes  T  long  to  modulate  the 
%  received  signals  from  each  antenna 

%Construct  a  PN  object  usnign  MATLAB  pn  seq  gen 
%  h  =  seqgen. pn(' Shift ' ,  0); 

%  %  Output  chipped  PN  bits 
%  set(h,  ' NumBitsOut ' ,  Nchips); 

%  for  i=l:N 

%  B(i,:)  =  generate(h); 

%  end 

%  B  =  Dmap(B,l);  %  map  0,1  to  +-1 

%  use  Walsh  Sequence  instead  of  pure  pn 

H  =  hadamard(Nchips) ; 
for  i=l:N 

B(i, :)=H(i+l, %  don't  want  1st  row 
end 

%  spread  B  by  SPC 

B  =  repmatC  B,SPC,1);  %  repeat  by  spread  factor 

B  =  reshape(B,N,SPC*Nchips);  %  like  chips  together 
B  =  B*pi/2;  %  pi/2  since  inside  exponential 

%test  SDMA  concept  by  elimnating  the  pn  sequence 
%B=zeros(N , SPC*Nchips) ; 

test=-30*pi/180; 

test2=45*pi/180; 

Vr=zeros(361,length(t)) ; 

Vr2=zeros(361,length(t)) ; 
er=zeros(361,length(t)) ; 
for  k=l : 361 ; 

an(k)=(k-l)-180; 

th(k)=an(k)*pi/180; 

%  calculate  the  actual  received  signal  for  each  angle  using  the  N  element 
%  receiver  and  sum  all  N  signals, 
for  i=l:N  %  create  combined  signal  at  angle  k 

Vr(k,:)=  Vr(k,:)+  exp(lj*Czeta+2*pi*Ci-l)*d*sin(test)  +B(i,:))); 
er(k,:)=  er(k,:)+  exp(lj*(  2*pi*Ci-l)*d*sin(th(k))+BCi , : ))) ; 
Vr2(k, :)=Vr2(k, : )+exp(lj*Czeta+2*pi*Ci-l)*d*sin(test2)+BCi , :))); 


end 

Vr(k, :)=Vr(k, :)+Vr2(k, 

RlCk)=sumCVr(k, : ) . *  er(k,:)); 

R2Ck)=sumCVr(k, : ) . *conj(er(k, :))); 
phasel(k, : )=unwrapCangle(Rl(k))) ; 
phase2(k , : )=unwrapCangle(R2(k))) ; 

end 

AF=sin(N*pi*d*(sin(th)-sin(test))+eps) ./CN*sinCpi*d*Csin(th)-sin(test)))+eps) ; 

AF2=sin(N*pi*d*(sin(th)-sin(test2))+eps) ./CN*sinCpi*d*Csin(th)-sinCtest2)))+eps) ; 

maxAF=maxCabs(AF)) ; 

maxAF2=max(abs(AF2)) ; 

maxR2=max(abs(R2)) ; 

maxRl=maxCabs(Rl)) ; 

figure; 
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plot(an,abs(R2)/maxR2, 'b' ,an,abs(AF+AF2), 'g') 

%  plot(an,abs(R2)/maxR2, 'b' ,an,abs(AF), 'g' , an, abs(AF2), 'r') 

xlabel(  \theta  ) 

ylabel(  Correlation  Magnitude') 

titlestring0  =  ['Uniform  Linear  Array  -  SDMA']; 

titlestring=['N=' ,num2str(N), '  d=\lambda/2, ' , '\theta  inc  =  '... 

,num2str(test*180/pi), ' ,  ,num2str(test2*180/pi)  ]; 
titlestring2=[' Nchips= ' ,num2str(Nchips), '  Samples/chip  =',  num2str(SPC)] ; 
title({titlestringO, titlestring, titlestring2}) 
axis([-90  90  0  1]) 
legend(  Rk' , 'Array  Factor  ); 
grid  on 

%  figure; 

%  titlestring=[' Correlation  Comparison  \theta_k  =' ,num2str(test*180/pi)] ; 
%  titlestring2=' comparison  of  instantaneous  and  averaged  phases'; 

%  plot(t,zeta,t,phase2(121, : )) 

%  xlabel(‘time  (t)') 

%  ylabel('phase') 

%  axis(  [0  t(length(t))  -2  2]) 

%  grid  on 

%  title({titlestring,titlestring2}); 

%  legendC ' phase  zeta(t) ',' average  phase') 


E.  SDMA  Technique  Applied  to  the  Random  Planar  Array 

This  code  was  used  to  generate  Figure  24  the  Signal  DOA  determination  for  a 
random  planar  array. 

%  Chris  Taylor  Modification  to  Dr.  Gross'  SDMA  sim 
%  to  allow  for  Random  Planar  array  vice  ULA 
%  Feb  08  Naval  Postgraduate  School,  Monterey  CA 

% 

%  Algorithm  from  Frank  B.  Gross  "Smart  Antennas  for  Wireless 
%  Communications",  c  2005,  Chapter  8,  p253-260 

% 

%  This  is  a  simulation  to  model  the  Elam  SDMA  technique 
%  Developed  by  Frank  B.  Gross,  November  11,  2004 
%  WALSH  FUNCTIONS  WORK  BETTER  THAN  PN  SEQUENCES 
%  Can  also  use  Gold  codes 
clear;  clc;  elf  reset; 

%  -  Define  Constants  - 


define  the  sample  baseband  modulation  initially  as  a  300MHz  phase  modulation  tone 


f=3E8; 

T=l/(4*f) ; 
Nchips=64; 

SPC  =  16; 

num_elements  =  30; 
phi_inc=30*pi/180; 


freq  of  baseband  modulation 
time  bandwidth  product  =  4 
Number  of  Chips  /  Time  seq 
Samples  Per  Chip 
number  of  antenna  elements 
incident  signal  of  interest 


lambda  =  3E8/f; 
k  =  2*pi/Lambda; 


%  c/f 

%  wavenumber 


tend=SPC*Nchips ;  %  time  snapshot  end 

t=[0:(tend-l)]*T/(tend-l);  %  discritization  of  signal 

zeta=sin(2*pi*f*t) ;  %  modulation  zeta  (the  signal  phase  sent) 
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Ntheta=91;  Nphi=361;  %  Set  up  angle  matrices 

%Ntheta=5;  Nphi=13; 

theta  =  linspace(0, pi/2, Ntheta); 

phi  =  linspace(-pi,pi,Nphi); 

DTheta=  theta(2)-theta(l) ;  DPhi  =  phi(2)-phi(l) ; 

[Phi, Theta]  =  meshgrid(phi , theta); 

CT=cos(Theta);  ST=sin(Theta) ;  CP=cos(Phi);  SP=sin(Phi); 

Rx=ST.*CP;  Ry=ST.*SP;  Rz=CT ;  %  unit  radial  vector  components 

%  -  generate  pseudo-random  chip  codes  - 

%  T  long  to  modulate  the  received  signals  from  each  antenna 
%  "codes"  here  =  "Beta"  nomenclature  from  text 
H  =  hadamard(Nchips);  %  use  Walsh  Sequence 
for  i=l : num_elements 

codes(i, :)=H(i+l, %  don't  want  1st  row 
end 

%  spread  codes  by  SPC 

codes  =  repmatC  codes , SPC , 1) ;  %  repeat  by  spread  factor 

codes  =  reshape(codes,num_elements,SPC*Nchips);  %  move  like  chips  together 

codes  =  codes*pi/2;  %  pi/2  since  inside  exponential 


Build  Random  Array 


array_size  =  5.0*lambda;  %  in  wavelengths  (lambda)  square  grid 

min_spacing  =  0.1*lambda;  %  l/10th  wavelength 

sepn_min  =0;  %  initialize  prior  to  1st  test 

xp  =  array_size  *  rand(l,num_elements); 

yp  =  array_size  *  rand(l,num_elements); 

% - Stats  of  Random  Array - 

while  (sepn_min<  min_spacing)  %  min  spacing  =  .1  wavelength 

rp  =  [xp;yp] ' ; 

sepn  =  pdist(rp);  %  dist  between  all  elements  in  Random  Array 
sepn_min  =  min(sepn);  %  find  minimum  separation 

zsepn  =  squareform(sepn) ;  %  convert  linear  sepn  array  to  square 

[Rmin  Cmin]  =  find(zsepn==sepn_min,l,  'first');  %returns  index 
fprintf(['Min  separation  is  %6.3f  wavelengths,'... 

'between  elements  %2.0f  and  %2.0f  \n '],... 
sepn_min,  Rmin,  Cmin) 

if  (sepn_min<  min_spacing)  %  then  Move  one  of  them 
xp(Rmin)=array_size  *  rand; 
yp(Rmin)=array_size  *  rand; 

end 

end 

fprintf(  'Element  %2.0f  Located  at  %8.3f,  %8.3f\n' ,Rmin,rp(Rmin, :)  ) 
fprintf(  'Element  %2.0f  Located  at  %8.3f  %8.3f\n' ,Cmin, rp(Cmin, :)  ) 

% - p-[0t  Random  Array - 

subplot(2,2,l) 

set(gcf,  ' PaperOrientation' ,  'portrait  ); 

plot(xp,yp, 'MarkerSize ' , 15, 'Marker ' LineWidth' ,2, . . . 

'LineStyle' , 'none' , 'Color' , [0  10]); 
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tstrg  =  'Sensor  Cluster,  Uniform  Excitation'; 

tstrg2=  [num2str(num_elements) , '  Elements  ]; 
title({tstrg, tstrg2} ,  ' Font Size' , 14) 
axis([0  array_size  0  array_size]) 
axis  square 

set(gca, 'XTick' ,linspaceC0,array_size,6)) 
set(gca, 'YTick' ,linspace(0,array_size,6)) 
xlabelC  x  -  wavelengths') 
ylabel(  y  -  wavelengths') 

grid  on; 

for  i  =  l:num_elements; 

text(xp(i) ,yp(i) ,  num2str(i),  ' FontSize' ,10) 
end 


% - plot  coord  system - 

subplot(2,2,2) 

coord_syst_diag(l , phi_inc*180/pi , 90 , array_size) 
plot3(xp-array_size/2,yp-array_size/2,zeros(length(xp)), . . . 

'Marker' , 'LineStyle' , 'none' , 'Color' , [0  1  0]); 
plot3(CP(l, :),SP(1, :),  zeros(size(CP(l, :))), 'c  ) 
titleC Coordinate  System  and  Scan  Plane'); 


%  Calc  Array  Factor  for  Planar  Random  Array 

% _ 


%  Half-wave  dipole  pattern  for  Iin=lA  (eps  handles  theta=0  indeterminent) 
E0  =  60*cos(pi*CT./2)./(ST+eps); 

10=1;  %  assuming  uniform  current  amplitudes 

Ip=  10  *  exp(- j*k*(cos(phi_inc)*xp+sin(phi_inc)*yp)) ; 

%  Summing  far-field  pattern 

AF=zeros(size(E0)) ; 
for  n=l : num_elements 
xpRx  =  xp(n)*Rx; 
ypRy  =  yp(n)*Ry ; 

AF  =  AF  +  Ip(n)*  exp( j*k*(xpRx+ypRy)) ; 

end 


Calc  Received  Signals  via  SDMA 


yr=zeros(Nphi , length(t)) ; 
ye=zeros(Nphi , length(t)) ; 

m=find(theta==pi/2) ;  %  want  row  for  xy  plane 

Rxm  =  Rx(m,:);  %  lxNphi  vector 

Rym  =  Ry(m, : ) ; 


for  p=l:Nphi; 

%  calculate  the  actual  received  signal  for  each  angle  using  the  N  element 
%  receiver  and  sum  all  N  signals, 
for  n=l:num_elements 
xpRx  =  xp(n)*Rxm(p) ; 
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ypRy  =  yp(n)*Rym(p); 

re  =  cos(phi_inc)*xp(n)+sin(phi_inc)*yp(n) ; 

yr(p, : )=yr(p, : )+exp( j*(zeta+  k*(xpRx+ypRy)  +codes(n, :))); 

ye(p, : )=ye(p, : )+exp(j*(  k*(re)  +codes(n, :))); 

end 

R(p)=sum(yr(p, :).*conj(ye(p, :))); 

end 

maxR=max(abs(R)) ; 


Calc  Gain, 


Side  Lobe  Levels,  HPBW,  etc. 


%  Horizontal  pattern  (the  x-y  plane)  and  3dB  beamwidth 

n=find(phi==0) ; 
m=find(theta==pi/2) ; 

F2  =  abs((E0.*AF)/(num_elements)).A2;  %  Normalized  power  pattern 
IntF2=DTheta*DPhi*sum(sum((F2) . *sin(Theta))) ; 

Gain=4*pi*F2/IntF2; 

Gmax=max(max(Gain)) ;  GmaxdB=10*logl0(Gmax) ; 

%  Convert  Correlation  Magnitude  to  Gain  Pattern 
RF2  =  abs((E0(m, :).*R)/(num_elements)).A2; 

GainR2  =  4*pi*RF2/IntF2;  %use  same  normalizing  factor 
GRdB  =  10*logl0(GainR2) ; 

%  Steinberg's  Random  element  Gain  pattern  Eqn  8.10 
F2_Rand_ave  =  F2.*(l-l/num_elements)+l/num_elements; 

Gain_Rand_ave  =  4*pi*F2_Rand_ave/IntF2; 
GRA_dB=10*logl0(Gain_Rand_ave(m, : )) ; 


G3dB=Gmax/2; 

GH=Gain(m, : ) ; 

mf=find(GH  >=  G3dB,l, 'first'); 
mp=find(GH  >=  G3dB,l, 'last'); 

BWH=180*Cphi(mp)-phiCmf))/pi ; 

GHdB=10*logl0(GH); 

AFmax  =  max(absCAF(m, : ))) ; 

% -  cq-lc  Side  Lobe  levels  - 

BW  =  imregionalmax(GHdB) ;  %MATLAB  rtn  that  finds  local  maxima->BW  array  "l"s 

A1  =  find(BW);  %index  pointers  to  local  Maxima 

[bb,ix]=sort(GHdB(find(BW)), 'descend');  %  sort  by  highest  peak->lowest 

SLL=  bb(l)-bb;  %  Side  Lobe  levels,  down  from  peak 

SLi=  Al(ix);  %  index  to  peaks 

SLa=  180*phi(SLi)/pi ;  %  Angle  where  peaks  occur 

fprintf( ' First  Sidelobe  Level:  %6.2f  below  max  at  phi=  %4.0f  \n' ,SLL(2),SLa(2)) 
fprintf( ' Second  Sidelobe  Level:  %6.2f  below  max  at  phi=  %4.0f  \n' ,SLL(3),SLa(3)) 
fprintf('Mean  Sidelobe  Level:  %6.2f  below  max  \n' ,mean(SLLC2:length(Al)))) 


% _ 

% 

%  Plot  Array  Factor  vs.  SDMA  Pattern 

% _ 
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%  subplot(2,2,3:4) 

%  plot(phi*180/pi , abs(R)/maxR, ' b' , phi*180/pi , abs(AF(m, : ))/AFmax, ' g' ) 

%  xlabel( '\phi ' ) 

%  ylabelC Correlation  Magnitude') 

%  tstring=[ ' Number  Elements=' ,num2str(num_elements), . . . 

%  '  Nchips=' ,num2str(Nchips)] ; 

%  tstring2=['\phi  inc  =  ' ,num2str(phi_inc*180/pi),  '  M=',  num2str(length(s))] ; 

%  title({tstring,tstring2} , 'FontSize' ,14) 

%  axis([-180  180  0  1]) 

%  legendC'Rk' , 'Array  Factor'); 

%  grid  on 


%. 


%  Plot  Gain  vs.  SDMA  Gain  Pattern 

% _ 


GHMax=max(GHdB) ; 

GRMax=max(GRdB) ; 

GRAMax=max(GRA_dB) ; 

Diff=GRMax-GHMax;  %  SDMA  only  calculated  in  one  planar  cut  vice  sphere,... 

%  Normalization  factor  is  adjusted  to  acct  for  diff. 

phid=phi*180/pi ; 
subplot(2,2,3:4) 

plot(phid,GRdB-Diff , 'b' ,phid,GHdB, ' g ' , 'LineWidth' ,1.5) 

%plot(phid,GHdB, 'b' ,phid,GRdB-Diff , 'g' ,phid,GRA_dB, 'r' , 'lineWidth' ,2) 

v=axis;  v(l)=-180;  v(2)=180;  if  v(3)  <  -60,  v(3)=-60;  end 

axis(v) 

tl  =  ['Horizontal  Plane:  G_{max}  =  ' ,num2str(GmaxdB, '%8.1f'), 'dB']; 
t2  =  ['HPBW=  ',  num2str(BWH, ' %8. If '), '  deg',... 

'  SI L= ' ,  num2str(S!t(2) ,  "%6. 2f  below  Peak')]; 
title({tl, t2} , ' FontSize' , 14) 
xlabelC  \phi  (deg) ' , ' FontSize ' , 14) 
ylabelC  Gain  Pattern  (dB)' ,' FontSize' ,14) 
legendC'SDMA  Correlation',  'Planar  Array  Gain'); 
grid  on 
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F.  Incident  Signal  Phase,  Walsh-Hadamard  Codes  &  Received  Signal 

This  code  used  to  generate  Figure  22  showing  the  received  signal  phase,  Walsh- 
Hadamard  Spreading  Codes  and  the  real  part  of  the  chipped  signal  r(t). 

%  This  is  a  simulation  to  model  the  Elam  SDMA  technique 
%  Developed  by  Frank  B.  Gross,  November  11,  2004 

%  define  the  sample  baseband  modulation  initially  as  a  1kHz  phase  modulation  tone 
%  WALSH  FUNCTIONS  WORK  BETTER  THAN  PN  SEQUENCES , Can  also  use  Gold  codes 

%  Modified  by  Chris  Taylor  Jan  2008  to  include  pn  seq  generation  or  Walsh 
%  codes,  samples  per  chip,  and  two  incident  beams 

%  This  short  code  shows  the  phase  of  incident  signal  zeta(t),  Walsh  codes 
%  and  the  Real  part  of  the  received  signal  r_n(t) 

clear;  clc; 

f=lE3; 

T=l/(4*f); 

Nchips=64; 

SPC  =4;  %  Samples  Per  Chip  (4  min) 

tend=SPC*Nchips ;  %  time  snapshot  end 

t=[0: (tend-l)]*T/(tend-l) ; 
zeta=sin(2*pi*f*t) ;  %  modulation  m 

%  plot  array  factor  for  linear  array 

N=ll;  %  number  of  receive  elements  to  give  10  degree  beamwidth 
d=.5;  %  spacing  in  wavelengths  between  adjacent  array  elements 

%  generate  pseudo-random  chip  codes  T  long  to  modulate  the 
%  received  signals  from  each  antenna 

%  use  Walsh  Sequence  instead  of  pure  pn 

H  =  hadamard(Nchips) ; 
for  i=l:N 

B(i, :)=H(i+l, :);  %  don't  want  1st  row 

end 

%  spread  B  by  SPC 

B  =  repmat(  B,SPC,1);  %  repeat  by  spread  factor 

B  =  reshape(B,N,SPC*Nchips);  %  like  chips  together 
B  =  B*pi/2;  %  pi/2  since  inside  exponential 

figure; 

subplot(3,l,l) 
plot(t,zeta, 'b  ) 
axisC  [0  T  -2  2]) 

xlabelC'time  (t)') 
ylabelC  Phase  -  \zeta(t)') 
grid  on 

subplot(3,l,2) 

plot(t,B(ll, :), 'r  )  %  use  arbitrary  WH  row 

axisC  [0  T  -2  2]) 

xlabelC’time  (t)') 

ylabelC 'Walsh-Hadamard  Code') 
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grid  on 


subplot(3,l,3) 

plot(t, real(exp(j*(zeta+B(ll, :)))), 'm') 

axisC  [0  T  -2  2]) 

xlabel(  time  (t)') 

ylabelC  Real  Part  of  r_n(t)  ) 

grid  on 
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